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Abstract 


This investigation is concerned with the dynamic synthesis of a heli- 
copter. The method of approach is a variation of the coniponent-niode syn- 
thesis in the sense that it regards the aircraft as an assemblage of inter- 
connected substructures. The equations of motion are derived in general 
form by means of the Lagrangian formulation in conjunction with an orderly 
kinematical procedure that takes into account the superposition of motion 
of various substructures, thus circumventing constraint problems. 
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1 . Introduction 

This investigation is concerned with the’ development of a mathematical 
model capable of simulating the vibrational characteristics of a helicopter 
in various flight regimes, such as hover and forward flight. The helicopter 
represents a very complicated structure consisting^of a given number of inter- 
connected substructures. The object is to produce a mathematical formulation 
which includes all the system dynamic characteristics and yet is not so cum- 
bersome as to defy analysis. 

The method of approach represents a variation of the component-mode 
synthesis (Refs. 1 and 2) in the sense that it regards the aircraft as an 
assemblage of interconnected substructures acting as a system. The various 
substructures identified are the airframe, the transmission shaft, the main 
rotor (consisting of the hub and the rotor blades), and the tail rotor (re- 
garded as a rigid fan). To ensure that the various substructures are acting 
as parts of the whole structure, an orderly kinematical procedure is de- 
veloped which takes into account automatically the superposition of motions. 
This procedure does away with the question of constraints. 

The approach based on the substructure concept has the advantage that 
it permits a large measure of versatility in the mathematical modeling of 
the substructures. For example, due to the complex configuration, the air- 
frame is best represented by a discrete model. On the other hand, the trans- 
mission shaft or a rotor blade can best be represented by a continuous model. 
In the final synthesis, each of the substructures is simulated by only a 
limited number of degrees of freedom. To this end, one expresses the dis- 
placements as a superposition of space-dependent modes multiplied by time- 
dependent generalized coordinates. For the space-dependent functions or 
vectors one must use rigid-body modes and deformation modes capable of des- 
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cribing the motion of the substructure with sufficient accuracy. This 
Rayleigh-Ritz type approach permits the simulation of the aircraft by a dis- 
crete system. 

To derive the system equations of motion, it is desirable to use an 
approach devoid of unnecessary complications. An approach fitting this 

V, 

description is the Lagrangian approach which is capable of producing the 
equations of motion of the system without the need of calculating constraint 
forces acting at points connecting various substructures. To this end, one 
must use a consistent kinematical representation to calculate the inertial 
velocity of every mass point of the aircraft. Such a representation necessi- 
tates the introduction of a number of reference frames associated with the 
various substructures. The motion of these reference frames can be regarded 
as representing the rigid body motion, and the motion of a point relative to 
the frames can be regarded as representing the elastic motion of the sub- 
structure. 

The Lagrangian approach requires the calculation of the kinetic energy, 
the potential energy, and the nonconservative virtual work for the elastic 
system. The potential energy is due to the structural elasticity and gravi- 
tational forces. On the other hand, the nonconservative virtual work is due 
mainly to aerodynamic forces. A relatively complicated procedure is necessary 
to transform the physical aerodynamic forces and torque into generalized 
forces . 

The equations of motion are nonlinear and their general solution is 
beyond the state of the art. Fortunately, general solutions are not necessary, 
and in fact may not even be very informative. Indeed, quite often the in- 
terest lies in solutions in the neighborhood of certain particular solutions, 
where the latter are generally known as trim solutions. Trim solutions of 
special interest are those associated with hover and forward flight. To 
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examine the motion in the neighborhood of trim solutions, it is necessary 
to derive the so-called variational equations, which are really perturbation 
equations about the trim. 

In deriving variational equations, one must make certain simplifying 
assunptions. These assumptions require a judgement- as to the expected or- 
der of magnitude of the various generalized coordinates and system parameters. 
Such a set of assumptions constitute a so-called ordering scheme and it 

applies primarily to the main rotor. 

The derivation of the variational equations implies an extremely large 
nunfcer of matrix multiplications and differentiations with respect to the 
state variables and time. In addition, it implies the selection from the 
multitude of terms only those which the ordering scheme deems essential. 

The magnitude of the task demands an efficient approach. To this end, a 
procedure for the derivation of the variational equations by means of com- 
puter manipulation is developed. 


2. Generalized Coordinates and Velocities 


We shall be concerned with the flight of a helicopter in the neighbor- 
hood of a given point on the earth's surface, so that the gravitational 
field can be regarded as uniform. It will prove convenient to introduce 
a set of inertial axes XYZ with the origin at a point 0 on the surface of 
the earth, so that axes X and Y are in the local horizontal plane and Z 
is aligned with the local vertical, and measure the motion relative to this 

inertial frame. | 

The mathematical formulation will be produced by regarding the heli- | 

copter as an assemblage of a given number of interconnected substructures. | 

This requires a set of kinematically consistent coordinates, in the sense 
I that the motion of one substructure must take into account the motion of 

another. A kinematically consistent formulation can be obtained by des- 
cribing the motion of the substructures independently and then imposing con- 
straints guaranteeing that points shared by two substructures undergo the 
same motion. It can also be obtained by means of an orderly kinematical 
I procedure taking into account automatically the superposition of motions. 

[ For example, if the airframe, the transmission shaft, and the rotor are iden 

! tified as different substructures, then the absolute motion of the rotor can 

I be regarded as a superposition of the motion of the airframe relative to the 

inertial space, the motion of the transmission shaft relative to the air- 
frame, and the motion of the rotor relative to the transmission shaft. This 
procedure eliminates the need of constraint equations and is the one used 
in this investigation. 

To describe the motion of the airframe, it will prove advantageous to 
Introduce a set of airframe body axes x^V/^z^ origin A at an arbi- 

trary point in the undeformed airframe, so that x. is along the forward 




direction, is along the vertical, and is normal to both so as to 
form a right-hand system. Then the motion of the airframe can be des- 
cribed in terms of the translation of the origin A, the rotation of the 
airframe body axes relative to the inertial axes XYZ, and the elas- 
tic motion of the airframe relative to The translation of A rela- 

tive to XYZ is given simply by the position vector Wq^ from 0 to A and 
the absolute velocity vector Wq^. They can be expressed in terms of the 

column matrices (Wq^) = [wo„«oAY’'oAZ^^ ‘*0a’ ' ^”0AX“0AY*0AZ^^’ 
pectively, where the meaning of the vector components is obvious. On the 

other hand, the rotation of axes x^Y^^z^ relative to XYZ is fully determined 
by the matrix of direction cosines between these two sets of axes and by 
the angular velocity vector of axes x^y^z^ relative to the inertial space. 
In defining the orientation of axes x^y^z^ relative to axes XYZ, it is con- 
venient to regard x^y^z^ as a triad originally coincident with axes XYZ and 
moving relative to these axes. Then, the orientation of axes x^y^z^ rela- 
tive to axes XYZ can be obtained by means of three rotations: about z^, 

about x^, and about y^ in that order (see Fig. 2). This permits us 
to write the coordinate transformation from one set of axes to the other in 
the compact matrix form 

"•a’ ■ tW'"* <” 

where 

(Ca) • [x^ . (R) . [X Y zf (2) 

are position vecv. '"s and 


[T,„] . [x,][x,][x^] 
is the transformation matrix, in which 


(3a) 
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— 1 
o 

0 “ 


fex^ 0 -SX ‘ 


■ cx^ sx^ 0‘ 

eg- 

0 


. [Xy] • 

y y 

0 1 0 

. [x^3 - 

-sx^ cx^ 0 


Lo -SX, 



-sxy 0 cxy. 


. 0 0 1. 


(3b) 


Note that Eqs. (3b) represent transformation matrices defining the rotations 
described above. In which sx^ - sin x^^. cx^ » cos etc. Also from Fig. 2 

It can be easily verified that the absolute angular velocity vector 
which defines the angular velocity of axes x^y^z^ relative to axes XYZ and 
which Is also equal to the relative angular velocity vector can be 
written In the matrix form 

{Ra) ' {o»^} = x^[Xy][x^]{e3) + Xj^[Xy]{e^) + yeg) (4) 

where 



are unit vectors written in matrix form. 

Equation (4) gives in terms of components along axes x^y^z^, 
whereas {Wq^} » C%x'^0Ay'^ 0AZ^^ expressed In terms of Inertial components. 
However, there are equations requiring that {Wq^^} and (n^) be expressed In 
terms of components along the same system of axes. It Is often convenient 
to write {Wq^} In terms of components along axes x^y^z^, which can be done 
by simply premultl plying (wq^) by the transformation matrix [T^q]* 

The elastic motion of the airframe requires first a definition of the 
airframe Itself. The major question Is whether the power train, consisting 
of the transmission assembly, the transmission shaft, and the control 
linkages, should be considered as part of the airframe or as a separate 
subsystem. Moreover, there Is the question of the tall rotor. In view of 
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the difficulties encountered in accounting properly for the flexibility of 
the transmission shaft in airframe free-free modes, and because we desire 
a formulation independent of the type of airframe modes used, the trans- 
mission assembly and shaft will be treated ?»s a separate elastic subsystem 
connecting the airframe to the main rotor hub. As such, we will ignore the 
mass of the transmission itself* and account for the transmission mount 
flexibility by introducing two torsional springs at the base of the shaft 
acting along orthogonal axes which are perpendicular to the shaft. The 
flexibility of the transmission shaft itself will be described in terms of 
two bending and one torsional displacement. The tail rotor will be re- 
garded as a rigid fan, not part of the airframe, so that the airframe is 
assumed to exclude the transmission, transmission shaft, control linkages, 
and tail rotor. Note also that the control linkages are not being con- 
sidered explicitly but only through their kinematical effects. 

Next, let us consider any arbitrary mass point in the airframe. The 
position of this point relative to the inertial space is given by 


' '^OA ^ 


rA " ^A 


( 6 ) 


where r^ is the radius vector from A to the point in question when the 
airframe is undeformed and u^ is the elastic displacement vector of that 
point. Recognizing that and u^ are measured relative to the moving 
axes x^y/^z^t the inertial velocity of the point is simply 


■ 'ioA * ^A-f^A * «A> * % 


where Wq^ and were defined earlier and 


■ ^OA - <^A " “a’-^A * ^A 

u^ is the elastic velocity vector 


♦This assumption may be discarded later if it is 'ound that the inertia of 
the transmission is significant. 
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of the mass point measured relative to Equation (7) can be written 

in the matrix form 


<«A> “ CTaO^‘“OA> - I^'a " “a^'^'a’ * ‘“A> 
where a skew symmetric matrix defined by 


(8) 


[’"a * “ 


"Az * \z 

L-f'-Ay * %> 


■<’‘A2 * “Az' 


r. + u, 
Ax Ax 


% * “Ay 

-<% * v> 


(9) 


in which r^^, r^^, r^^ are the components of and u^^ are the 

components of u^. 

We note that the vectors {Wq^} and {n^} represent the translation of 
the origin A of axes x^y/^z^ aod the rotation of these axes. They can be 
interpreted as the "rigid-body modes" of the airframe. On the other hand, 

{u^} represents the elastic displacement vector relative to the airframe 
axes x^y^z^. As customary in the analysis of complex structures, we assume 
that the elastic displacements of the airframe can be represented as a 
linear combination of space dependent functions multiplying time dependent 
coordinates, where the first are referred to as "airframe modes". Th:re 
are two main possibilities. One is to use "free-free modes", which can be 
obtained by regarding the airframe as being clamped at point A, where A is 
taken to coincide with the center of mass of the undeformed airframe. Another 
possibility is to use "cantilever modes"* which can be obtained by regarding 
the airframe as being fixed at the base of the transmission. The implica- 
tions of these possibilities will be examined in the next section. 

Equation (8) defines the motion of any mass point of the airframe. One 
point of particular interest is that corresponding to the lower end of the 
transmission shaft, namely, the end attached to the transmission and trans- 
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trission mount. We shall denote this point by S and attach a set of axes 
^ES^ES^ES airframe at this point such that coincides with the 

shaft axis before any deformation takes place and x^^ and y^^ are attached 
to the airframe and are normal to z^^ (Fig. 3). Another set of axes, the 
shaft axes x^y^z^, have the origin at the same point S and are obtained from 
axes x^^y^gZ^^ by means of three rotations. Considering axes x^y^z^ to be 
initially coincident with axes x^^y^jZ^g* these rotations are i|)^ about x^, 

about y-, and ^ about z^. in that order, where if» , i() are small angles permit- 

y j z ^ X y 

ted by the transmission mount flexibility and = n is the constant angu- 
lar velocity imparted to the shaft by the engine. The position of point 
S relative to point A is given by the vector r^^ + u^^, where r^^ is the 
rcdi'js vector from A to S when the airframe is undeformed and u^^ is the 
elastic displacement vector of point S relative to axes x^y^z^. The trans- 
lational velocity vector w^^ of point S is obtained in matrix form by 
simply introducing the coordinates of point S in Eq. (8). On the other hand, 
the angular velocity vector of the frame x^y^z^ is simply 
where ^he elastic angular velocity vector of axes x^^y^gZ^^ relative 

to axes x^y^z^ and u)^ is the angular velocity vector of axes x^y^Zg rela- 
tive to axes x^^y^gZ^g. 

The mot. on of point S is in terms of components along axes x^y^z^. In 
particular, the translational velocity of S is 


("as’ " '■’^A0^("0a’ ■ ^■'as * “as^’^a’ * (“as’ 


no) 


where [r^^ + u^^j] can be obtained from Eq. (9) by replacing r^^ + (x^, 

yA» Zfl) by r^ 5 ^ + Ufl^^ (Xflc. y/ic' Zac)» etc., and where { 11 ^^} is a vector 


'ASx '^'AS’ -^AS’ ‘AS 


AS- 


having the components '^/\sy’ “aSz* interest lies in working 

with components along Xg, y^, z^, then these can be obtained by premultiplying 
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the components along by the transformation matrix [T^^] = 

representing the matrix of direction cosines between these two sets of axes, 
where [i^f5]C<iQ5] represents the matrix of direction cosines between axes 

^ES^ES^ES matrix of direction cosines between axes 

x^sYes^ES ^S^S^S* represents the matrix of direction 

cosines between axes x^y^z^ and x^^y^^z^^ due to the geometrical configuration 
when the airframe is undeformed, whereas is due to the airframe elas- 

tic deformations. To obtain the latter matrix, we assume that the orientation 
of axes x^^y^jZ^j is defined by three infinitesimal rotations given by the 
curl of the elastic displacement vector Therefore, we can write 


in the form 



u 


where u*^^, uj^^ are the components of the vector 

it is understood that the partial derivatives are to be evaluated at the 

point S. Furthermore, the matrix [t^] can be shown to have the explicit form 

= [<i'2]C'l'y]Ux] ( 12 a) 


in which 


[♦,] > 

1 0 
0 

0 

«• 

1 — 1 
V: 

11 

C|<y 0 -S4-y 

0 1 0 

. C’1’2^ = 

-Sil*^ 

0 

0 


^ -s*x 



1 V> 
-€• 
V: 

o 

o 


1 

o 

0 1 
(1 
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The angular velocity vector of the frame relative to 

is recognized simply as the curl of the elastic velocity vector u^^ which 
can be shown to have the vector form 

-ES ~ ^ ^AS 


where v is the well known del operator. Equation (13) can be written in 

•>s 

the matrix form 


in which 


(U) 


[V] = 


3Z 


ES 


ay, 


ES 


9Z 


ES 


3X 


ES 


(15) 


[ ay£s ^^ES J 

is the skew symmetric matrix operator, which is the matrix counterpart of 
the vector operator. On the other hand, the angular velocity of axes 
^S^S^S ^ES^ES^ES shown to be 


{(Og) = 

where {e.} (i = 1,2,3) are given by Eqs. (5). Hence, the total angular 
velocity of axes x^y^z^ has the form 

{fig} = ^ ^“S^ 


Next, we wish to define the motion of the shaft. To this end, let 
us consider an arbitrary point on the shaft. The position of this point 
relative to axes x^y^z^ is given by the radius vector 1[;5 + where r^ is 
the position of the point in question when the shaft is undeformed and u^ 
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is the elastic displacement of that point relative to axes Denoting 

by Ug the elastic velocity of the point relative to axes x^y^z^, the iner- 
tial velocity of this point can be written as 


“s ' V * ?S ■“ <rs * Ss> ^ “s ° “as ■ <r$ * “s> ” ?S * “s 


08 ) 


But is given in terms of components along axes-oc^y^z^, whereas the re- 
maining two terms are in terms of components along axes x^y^z^. To express 
all terms in Eq. (10) in terms of components along axes x^y^z^, we premultiply 
{w^^} by the transformation matrix [T^^], so that Eq. (18) can be written 
in the matrix form 


‘“s' ° tV‘“As' ■ '^''s * “s^'^s' * ‘“s’ 


(19) 


To describe the elastic motion, it is convenient to regard the shaft as 
a one-dimensional member. Note that this is in direct contrast with the 
airframe which was regarded as a three-dimensional structure. We can write 
the position of any point on the deformed shaft relative to axes x^y^z^ in 
the matrix form (Fig. 4) 


‘'"s “s’ 



(20) 


where u^, v^ are elastic displacements parallel to axes Xg, y^ respectively, 
of a point originally on the axis z^ and at a distance z^ from S. Axial 
elastic displacements have been assumed to be smaller than u^ or v^ and 
have been ignored. The shaft also undergoes torsion, but this does not 
affect {Wg}, which, in view of Eq. (20), is the velocity of a typical point 
on the shaft, originally on the axis Zg. The torsion does, however, affect 
the orientation of a set of axes attached to the shaft at any point z^ 
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which were originally parallel with axes and are moving with the shaft 

during deformation. Denoting by [T^l the transformation matrix from the 
undeformed shaft axes to the deformed axes, we can write 






'X. 


( 21 ) 


in which (|)^ is the elastic torsional displacement about axis and where 
primes denote differentiations with respect to the spatial variable z^. 

To define the motion of the rotor, let us assume that the upper end 
of the shaft coincides with the geometric center H of the hub. Then let 
us introduce the hub axes with the origin at H and with axis z^ 

along the rotor spin axis; axes x^ and y^ are attached to the hub and are 
normal to z^. Note that axes x^ and y^ are parallel to axes xg and yg 
when the shaft is undeformed. The relation between the direction of axes 

general form 


(r„) = [T^jKrj) (22) 

where the matrix [T^^] of direction cosines is simply the matrix [T^] given 
by Eq. (21) evaluated at z^ = L^, where is the length of the shaft. More- 
over, the angular velocity of the hub axes x^y^^z^ is 


where 


(23) 


*“Eh’ “ ijO-j.t)]^ (24) 

is the angular velocity of axes x^yuZu relative to axes x y z due to the 
elastic angular motions at the upper end of the shaft. The absolute velocity 
of point H is simply Eq. (19) evaluated at z^ = and can be written in 
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(25) 




the symbolic matrix form 

(«SH> = C"sa]‘«AS> - " <“SH> <''> 

The tail rotor is assumed to be a rigid fan spinning at the angular 
velocity relative to the airframe. Using the analogy with the transmission 
shaft, let x^jy^jZgj be a set of axes attached to the frame at the point T, 
coinciding with the center of the fan (Fig. 5 ) snd let XjyjZ^ be a set of 
axes rotating relative to translational velocity of 

T is 

<“AT> ' t^A03'"0A> - f'-AT " “aT^'“a> " ‘“AT> 
where all the quantities are as in Eq. (10) except that the coordinates of 
point T replace those of point S. Similarly, the angular velocity of axes 

^ET^ET^ET ’"®^3tive to axes x^y^z^ is 
= [vKUejKu^j}) 

and the angular velocity of the tail rotor is 

{fij} = ^ 

where there are no counterparts of and i|»y for the tail rotor and the 
spin axis is taken as coincident with axis Zj. 

The main rotor is assumed to have M identical articulated blades 
where the blades are assumed to have the flap-lag-pitch configuration (Fig. 


The subsequent derivations are concerned with a typical blade i (i = 1, 
2,...,M). Ordinarily, we would identify every quantity pertaining to the 
blade by the subscript i. However, with the tacit understanding that the 
subscript is implied throughout, considerable simplification of notation can 
be achieved by omitting the subscript during the derivation stage and rein- 
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troducing It when it becomes necessary. Hence, let us consider a typical 

blade and assume that the flap hinge Is at point F at a distance f>^om H, 

'Hr 

the lag hinge is at point L at a distance rp^ from F, and the pitch hinge 
is at point B at a distance from L. First, let us introduce the set of 
axes XpypZp obtained through a rotation og about axis and a rotation -e 
about axis y^, where ag is known as the azimuth angle and B is known as the 
flapping angle (Fig. 6a). Note that Og is constant for every blade. The 
relation between axes XpypZp and Xp^y^^z^^ is given by 

(29) 




where 


[Tph] = [B][oig] 

is a transformation matrix, in which 



Cctg 

sag 

o' 


C6 

0 

SB" 

[aj] = 

-sag 

Cag 

0 

. [6] = 

0 

1 

0 


0 

0 

1 


-SB 

0 

CB 


(30) 


(31) 


Although it is more convenient to define the components of the radius 
vector rpip in terms of components along an intermediate set of axes, for 
consistency of notation, we express r^p in terms of components along axes 
write 


The inertial velocity of point F has the matrix form 


(32) 


{Whf> = ‘ 

wh. > [r^^p] is the skew symmetric matrix corresponding to irp^p). Also from 
Fig. 6a, it can be verified that the angular velocity of axes XpypZp has 
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the form 


{fip} = [Tp^]{J2|^) * 6 ( 62 } (34) 

Following the same pattern, we conclude from Fig. 6b that the relation 
between the lag axes and the flap axes XpypZp is 




(35) 


where 


[Tlf] = 


Ca So. 0 

-Sa Ca 0 

0 0 1 

in which a is known as the lag angle, 
the lag hinge is 


(36) 


Moreover, the inertial velocity of 


{Wp^} = [Tp^]{w^p} - [rpj_]{fip} (37) 

in which [i^pl] is the matrix obtained from {>"p|_} = Lp{e.|}, where the latter 
is the vector from F to L whose magnitude is Lp. The angular velocity of 
axes X|^y|^Z|^ is 

Finally, we wish to define the blade axes XgygZg such that Xg is along 
the axial direction of the blade, yg is the plane of the blade and in the 
direction of the leading edge, and Zg is normal to the blade (Fig. 6c). The 
relation between axes XgygZg and X|^yj^Z|^ is 


(39) 


where 


(rj) = [TB|_](r,_) 


1 0 0 

0 ce se 

0 -se ce 

In which e is known as the pitch angle. 




(40) 

The Inertial velocity of point B Is 


{Wlb^ = 


where C'^Lg] "^s the matrix obtained from { 
vector from L to B whose magnitude Is L|^. 
axes XjygZj 1s 


(41) 

*^LB^ = where {r^^g} Is the 

The angular velocity of the blade 


{J2b> “ [TBL]{nL> + 0(ei} (42) 

At this point, we note the Implicit assumption that the hub and Its 

connecting links are rigid, with the connecting links being considered as 

one-dimensional bodies. In addition, we shall assume that the connecting 

links are small compared to the airframe and the rotor blades, so that they 

can be regarded as massless. Hence, there Is no direct contribution from 

these links to the system energy. In contrast, the hub can be much more 

massive, particularly for large helicopters, and hence will contribute to 

the system energy as a rigid body. 

The rotation angles og, B, a, and 0 described above deserve further 

discussion. In particular, the angle Ug Is a constant design parameter. On 

the other hand, the angles 6 and a representing flap and lag, respectively, 

are time-dependent generalized coordinates. The angle 0 Is the sum of three 

parts, namely, the collective pitch 0 , the cyclic pitch 0.^, and the 

pitch change 0 due to coupling, where 0^. and 0_„ are given functions of 
cp CO cy 

time and 0 Is a function of 6 and/or a and provides the flap and/or lag 
cp 

kinematical coupling, respectively. Hence, 0 does not Introduce additional 
degrees of freedom. 

We now wish to define the motion of a typical rotor blade. By analogy 
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WW r:-‘ 




with the airframe and the transmission shaft, we can write the absolute 
velocity of an arbitrary point on the rotor blade in the general vector 
form 


"b ' «LB * 8b ■< (rB * “b> * 8b ' «LB ■ '^B * 8b> '< “b * 8b 
where ro is the radius vector from B to the point in question when the blade 

-D 

is undeformed, Ug is the displacement vector of the point arid, Ug is the 
elastic velocity vector of the point relative to XgygZg. We note that w^^g 
is given in terms of components along axes x^y^z^ whereas the remaining two 
terms are in terms of components along x.^ygZg. Following the established 
pattern, we can write the velocity vector Wg in terms of components along 
axes XgygZg in the matrix form 


'«b> = [Tbl^'“lb> ' ^'•b * “b’’“b> * ‘“b> 

where [rg + Ug] is the skew symmetric matrix associated with rg + Ug. 

It remains for us to define the elastic motion of the blade relative 
to axes XgygZg. To this end, we can write the position of any point in the 
deformed blade relative to axes XgygZg in the matrix form (Fig. 7; see also 
Ref. 3) 



fx4u] 

/-x{(t»^ + ♦)•' 

(rB * Uj) = . 


■ ^ 


1 W i 

1 c 


where u, v, w are elastic displacements along axes Xg, yg, Zg respectively, 

e is the distance at the blade root between the pitch axis Xq and the 
0 ° 

elastic axis, measured in the Xgyg plane and positive in the direction of 
the leading edge, + <!»)' represents an axial position where x = x(n,c) 
is the warp function (note: x{0,0) * 0), and n, c are cross sectional coordi- 
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i 


^30*m 


nates along local principal axes at a distance x from B. in which n is along 
the chord pointing toward the leading edge and c is normal to the blade 
mid-plane. Moreover, [Tpg] is the transformation matrix from the undeformed 
blade coordinates to the deformed coordinates. 

When deriving the rotor blade kinetic and potential energy, it becomes 
necessary to omit certain small terms arising from'the blade elastic defor- 
mation in order to keep from unduly complicating our equations. These 
terms are generally known as "higher-order" terms. At this point, we wish to 
introduce the idea of an ordering scheme which will enable us to compare the 
relative magnitudes of terms in a systematic manner, thus permitting us to 
determine which terms to retain and which terms to omit from the formulation. 
To this end, let us define the nondimensional parameter c, which is taken 
to be the same throughout this report. The approximate magnitude of e is 
taken to be less than one tenth and a term of the same order of magnitude 
as c is said to be of order e and denoted 0(e). We can compare the relative 
magnitudes of terms in an equation by first nondimensional izing that equa- 
tion and then comparing the resulting nondimensional terms directly with 
the parameter e. Hence, to compare terms in Eq. (45), we divide Eq. (45) 
by the length of a rotor blade Lg and compare each resulting term to e. As 
in Ref. 3, we shall adopt the following ordering scheme for the blades 


^■O(e^) p-O(e) ^■■0(£) ?-= 0 (£) 

*•0 Lb Ln 


= 0 (e ) 4 , a 0 (e) v' « 0 (c) w' = 0 (c) 

•■B 

l^=0(c) ^^-=0(c) ♦Q^Od) 


0(c) 

■-B 

r=0(l) (46) 

‘'B 
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Note that primes denote differentiations with respect to the spatial variable 
X, and (|iq = ♦q(x), (|i = 4i(x,t) are blade pretwist and blade torsional dis- 
placement about the deformed elastic axis, respectively, where the pretwist 
is considered as a static torsional displacement. The angles ♦q and ^ are 
defined here differently than in Ref. 3 where the sum 4 >q + was taken to be 
the "total torsional bending" about an axis parallel to axis Xq. We shall 

D 

see in the next section that it is necessary to obtain the vector {r_ + Uq} 
given by Eq. (45) to 0(€'*). 

In view of Eqs. (46), the largest terms in the vector multiplying [Tpg]^ 
in Eq. (45) are 0(e) terms so that it is necessary to express [T^g] in terms 
of the elastic displacements to O(e^) only. Hence, the matrix [Tpg] can be 
written in the form (see Ref. 3) 





-[ v ' cos (^ Q +*)+ w ' sin (* Qn )] 

[ v ' sin ( 4 Q + t )- w ' cos (* Q +*)] 


(1 

(1 



v ' 

C 0 S (( J . Q +(|.+ v ' w ' ) 



) sin ( 4 > Q +( f .+ v ' w ' ) 


w' 

. w'2 

(1 — ^)sin(<ig+(|i) 
w'2 

(1 — 2“)cos((fiQ+^) 


(47) 


in which it is assumed that the pretwist and torsional displacement are both 
taken to be zero at the blade root. It is understood that consistent order 
series expansions are used for sine and cosine terms. Introducing Eq. (47) 
into Eq. (45) and expressing n and c in terms of components along yg and Zg, 
we obtain 




x+u-x((|.Q+^) '-(v'+w'^)(ncos*Q-csin^Q) - (w'-v>)(nsin«j>g+i;cosi)Q)j 

v'2 / 

eQ + v+(l — p )( ncos (|> Q - csin 4 > Q ) - ( i *+ v ' w ‘ )( nsin ((. Q + ccosiiQ ) 


..2 


w+(l - ^)(nsin(j>^+tcos*«) + ♦(ncos^n-^sintn) 


j 


( 48 ) 
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where O(c^) terms and smaller have been neglected. Equations (48) contain 
nonlinear terms which we desire to retain. This is in contrast to the air- 
frame and shaft, whose deformations are described strictly in terms of 
linear elasticity. 

Finally, it will prove convenient to obtain an expression for the angu- 
lar velocity of the deformed blade axes Xj^y^Zp at any point on the elastic 
axis a distance x from the blade root, where y^ and Zp are local principal 
axes of the cross section with origin on the elastic axis, in which y^ is 
along the chord pointing toward the leading edge and Zp is normal to the 
blade mid-plane, and where Xp is tangent to the deformed elastic axis posi- 
tive toward the free blade end. Note that n and c are measured along axes 
yp and Zp, respectively. An expression for the angular velocity of axes 

XgygZg can be deduced from the bending curvature 
expressions given in Ref. 3 where care must be exercized in distinguishing 
between the torsional displacement as defined herein and the elastic kine- 
matical pitch angle of Ref. 3 (see Ref. 4). Such an expression can be 
written in the form 


■ * I * 

w V ♦♦ 


=<(v'-vV^ - I v'w'^)sin^*-{w'-v'v'w'-w'w'^ - iwV^ 


)cos4>. 


(v'-v'v'^ - ^ v'w'^)cos*Q+(w'-vVw'-w'w'^ - wV^)sin^P 

+ v'(kos*q - j ♦^sin^pl+w'C^sinc^P + ^ ♦^cos^p) . (49) 

- v'(^sin^Q + j ♦^cos^q)+w'(*cos^q - ^ i>^sin4>p) 

where O(e^) terms and smaller have been neglected. Hence, the angular 
velocity of axes XpypZp is 


21 


REPRODUCBiLrry of the 

ORIGINAL PAGE IS POOR 



l«o) - HobKi'bI + (SO* 

where [T^g] is simply Eq. (47) evaluated at the point 0 in question. 

A summary of some Important equations in this section can be found 
in Appendix A. 


I 


2: 


3. The Kinetic and Potential Energy 


■ 


To derive Lagrange's equations of motion it is necessary to produce 
expressions for the kinetic energy, the potential energy, and the noncon- 
servative virtual work. In this section we shall consider the kinetic 
and potential energy for the system and in the next section we shall present 
the virtual work due to aerodynamic forces. 

The kinetic energy of the airframe can be written in the general form 



twAldmA 


(51) 


where mA is the mass of the airframe. Introducing Eq. (8) into Eq. (51), 
we obtain 




m. 


{UAldm^ 


+ twQA> [Taq] ^^AO^ ’ ^^A^ ^^A^ 


where 


<“a> 

[0,] 

<"a> 




[rA * 

(u^}dniA 

[■•a * “a’^[''a ^ “a^‘‘"’a 

tr-A * 


(52) 

(53a) 

(53b) 

(53c) 

(53d) 

(53e) 
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will be referred to as "mass integrals", where m^ is the total mass of the 
airframe, [J^] is the matrix of moments of inertia of the airframe in de- 
formed configuration about axes and {H^} is an angular momentum 

vector of the airframe due to elastic velocities. Note that Eq. (52) takes 
into account the fact that [T^q] is an orthonormal matrix. 

Next, let us consider the shaft kinetic energy" and write it in the 
general form 


T = - 
S 2 


m. 






where {w^} is given by Eq. (19). Introducing Eq. (19) into Eq. (54), we 
obtain 




iTtt tTi 


n I 1 T r 


where 


m^ = 


S Inr ""s 

s 


[fs] = 


m. 


[rj + Ujldirij 


(56a) 

(56b) 


(Us) = 


m. 


{Usldms 


[Jj] = 1 [fj + + Ujldm,. 


(56c) 

(56d) 


(Hj) = 


[r + u ]^{i )dm 

"s 


(56e) 


are certain mass integrals and is an orthonormal transformation. Equa 
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tion (55) looks deceivingly simple and is in fact quite complicated. For 
example, the vectors depend not only on the rigid body mo- 

tions of the airframe but on the airframe elastic motion as well. Moreover, 
[T$a] involves the elastic rotation of the airframe at S and the angles 
between axes ^nd which vary with time. In particular 

the translational velocity and angular velocity expressions, Eqs. (10) and 
(17) respectively, must be substituted into Eq. (55). Later in this sec- 
tion, we will simplify Eq. (55) by ignoring certain higher order terms. 

Assuming that the point T coincides with the mass center of the rigid 
tail rotor, the kinetic energy of the tail rotor can be written as 


Tj = -^ i (57) 

where m^ is the mass and [J^] is the matrix of moments of inertia about 
axes Xyy^z^P of the tail rotor. Similarly, assuming that the point H coin- 
cides with the mass center of the rigid hub, we can write the hub kinetic 
energy as 






(58) 


where m^ is the mass and [J^] is the matrix of moments of inertia about 
axes x^y^z^ of the hub. Although Eqs. (57) and (58) are considerably simpler 
than Eq. (55), the velocities (Waj} and {fjy} for the tail rotor and 
and for the hub still depend on the motion of the airframe. In addi- 
tion, {Wg^} and also include the motion of the transmission shaft. 

Finally, we can write the kinetic energy for each rotor blide in the 
general form 


T = — 

'b 2 


{Wg} {WsldmB 


(59) 
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where is the mass of the blade. Note that the subscript i, identifying 
a particular blade, has been ignored temporarily. Introducing Eq. (44) 
into Eq. (59), we obtain 


^ 2 "'b 


^ nin 


+ (w^^b)^[Tb|^]^(Uj) - (W|b)^[T„|] [rsDlnB' ■ 


LB' '•'BL-' '•'B" B‘ 


B B 


(60) 


where 


[rB] 

[Jb] 

(Hb) 




FTlp 


nir 


■I. 

■I. 


[••b ^ “b^‘^ 

<“b'*b 

[■•b ^ “b^^[''b “b^‘*"b 

^>■8 * " b ^^‘“ b >*^ 


(61a) 

(61b) 

(61c) 

(61d) 

(61e) 


The various quantities in Eqs. (60) and (61) are analogous to those per- 
taining to the airframe and the shaft. Of all the component kinetic energies, 
Eq. (60) is the most complex. In general, the velocity and angular 

velocity {^g} take into account the rigid body translation and rotation of 
the airframe, the airframe elastic motion, the rigid body rotation of the 
transmission shaft relative to the airframe, the transmission shaft elastic 
motion, the rigid body rotation of the flap link relative to the hub, and the 
rigid body rotation of the lag link relative to the flap link. A method 
of simplifying these coupling effects will be given later in this report. 

The blade ordering scheme of Sec. 2 will be used in this section to obtain 
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0(g ) expansions of Eqs. (61). 

Equations (52), (55), (57), (58), and (60) give the kinetic energy for 
each component. The system kinetic energy is formed by simply summing up 
these equations, where it is now necessary to introduce the subscript i in 
Eq. (60) and sum over the number of blades. This yields 

Substituting Eqs. (52), (55), (57), (58), and (60) into Eq. (62) and grouping 
together terms of the same type, we can write 
6 

T = I T. (63) 

j=l 


in which 

Ti = ^ rri;^ * I "’t^'^AT^'*'^'^AT^ \ 

1 • T • 1 ^ . T . 

+ 2 niH {W5^} {wj^} + 2 ^'^LB^i^\B^i 

■ ^\B^i^^'’^BL^i^’"B^i^“B^• 


+ 2 ’’’ 2 ^^B^i^'^B^i^^B^ 


^6 = 1 f ^ 1 1 + i J 

■'m. ■'m^ 3 3 


(64a) 


(64b) 

(64c) 


(64d) 

(64e) 


m, 


Bi 


(uBl^iuBlidmBi 


(64f) 
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Equations (64) possess all of the complexities mentioned earlier. A direct 
approach to evaluating the system kinetic energy explicitly is to substitute 
explicit expressions for Eqs. (53), (56), and (61) as well as explicit 
expressions for the velocities and angular velocities given by Eqs. (A28)- 
(A32) and (A12)-(A16) into Eqs. (64). The extreme tediousness involved 
forces us to seek explicit forms for some terms at 'one time and for other 
terms at a different time. This is motivated by the different treatments 
necessary for the two different types of quantities, those which depend on 
the spatial position and time and those which depend on time alone. One 
possibility is to substitute explicit expressions for Eqs. (53), (56), and 
(61) into Eqs. (64) while retaining all other quantities in implicit form, 
i.e., to write the components of the velocity and angular velocity vectors 
of each set of axes implicitly. This approach is convenient when the elas- 
tic displacements and velocities are of interest. It is the approach of 
this section. Another possibility is to substitute Eqs. (A28)-(A32) and 
(A12)-(A16) into Eqs. (64) while retaining the matrix quantities defined by 
Eqs. (53), (56), and (61) implicitly. Such an approach is more suitable 
when the discrete coordinates are being considered and the approach will be 
discussed in more detail later in this report. 

Let us now examine the airframe kinetic energy, Eq. (52), more closely. 
Assuming that the angular velocity vector and the elastic displacement 
vector sufficiently small that their product can be ignored, we can 

make the following simplifications in Eqs. (53b, d, and e) 


"a 

t^A^ ' I, 


(65a) 

(65b) 


(Ha> ' I 

™A 


(65c) 
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if. 

f: 

f: « . 

y-j 

f 

I Note that this is equivalent to ignoring the cross product x u^ in Eq. (7). 

I Equations (65) have simple physical meaning. Equation (65a) is the skew 

symmetric matrix associated with the vector from A to the airframe mass 
center in undeformed configuration, Eq. (65b) is the inertia matrix of the 
undeformed airframe about axes x^y^z^, and Eq. (65c) is an angular momentum 
vector due to the airframe elastic velocities. 

Thus far, we have been concerned with the airframe elastic displacements 
and velocities in a general way only. Now we wish to explore the possibility 
of representing these displacements in terms of airframe modes. Recognizing 
that these modes are three-dimensional, we can represent them by 3 x 1 column 
matrices (x/^» z^)) and write 

P 

h- "a> '»Al(*A> h- ^A»"Ai<‘> 

where n^^-(t) are generalized coordinates associated with these modes. Note 
that the mode gives the three displacement components at every mass 

point of the airframe. Introducing the 3 x P airframe modal matrix 

(67) 

and the airframe elastic generalized coordinate vector 
ir\p) = '’A2 • • * 

Eq. (66) can be written in the compact form 
tUA> = 

Substituting Eq. (69) into Eqs. (53) and taking into account Eqs. (65) we 
can write the airframe kinetic energy as 

TiFPRODUcmiLrry of thb 
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5 ^: 

I 

! 


Ta = 2 '"a * 2 ^^A^ ■*■ 2 ^^A^ ^^A^^^A^ 

+ ^Wqa)^[Tao]^[*A^^'^A^ ■ 


- <"a>'C‘a]'‘«a> 


(70) 

which we have introduced the definitions 



[M^] = [»^]^[»n]dm^ 


(7U) 

“A 



* 

[»^] =■ 


(71b) 

'""a 



t’A^ “ 1 t''A^^[*A^''"’A 


(71c) 

"■a 




Substantial simplification can be achieved if the point A is chosen to 
coincide with the center of mass of the undeformed airframe. Then, if the 
elastic modes are orthogonal to the rigid body modes, we can write 

° ^*A^*"a * f'^A^^'-*A^'*"A ° 

\ ™A "a 

SO that 

[rA] = [^a] = Ha^ = [0] (73) 

The requirement that ^aYaZa principal axes of the airframe does not 
yield significant savings, so that it need not be made. Inserting Eqs. (73) 
into Eq. (70), the kinetic energy reduces to 


where [Ja] is now the matrix of moments of inertia of the airframe in unde- 
formed configuration about axes ^a^a^A ^^A^ diagonal matrix of 
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"general! zed masses" . 

The above development is predicated on the availability of free-free 
airframe modes, which cannot be taken for granted. Moreover, it may not be 
the most advantageous way to describe the airframe motion. Hence, we wish 
to consider the possibility of using a different type of airframe modes. 

As mentioned earlier, an alternative to the use of airframe free-free 
modes is the use of cantilever modes, obtained by regarding the airframe 
as being fixed at the transmission base. We recall that this point was 
denoted by S in the preceding section. This alternative has the advantage 
that it eliminates the need of considering the intermediate axes Xr yr Zr- , 

* O " O " o 

and indeed one can assume that x^y^z^ take their place. On the other hand, 
the cantilever modes are not orthogonal to the rigid body modes, so that 
additional terms appear in the kinetic energy which do not appear in Eq. 
(74). 

In view of the above, let us assume that axes have the origin 

at point S and that axis z^ coincides with the direction of the shaft axis 
in undeformed state. Axes x^ and are attached to the airframe and are 
normal to z^. Then Eqs. (65a), (65b), and (71) remain valid except that 
their interpretation is different. Of course, the difference comes from the 
fact that the matrix [<t^] now represents cantilever modes and not free-free 
modes. Likewise, the vector {r^} is now measured from point S, which af- 
fects [<^^]» Cl^]» and [J^] . In particular, [J^] is now the inertia matrix 
of the airframe in undeformed configuration about axes x^y^z^ with the 
origin at S. Clearly, simplifications (72) are no longer possible and the 
kinetic energy remains in the form (70). 

In contrast to the airframe where all components of are small 
and vary with time, the third component of the shaft angular velocity vector 

contains a large constant part. Giving special consideration to this 
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angular velocity component, in the sense that quantities containing this 
component may not necessarily be small, we substitute Eq. (20) into Eqs. 
(56), integrate over the length of the shaft, and obtain 



(75a) 


(75b) 


[Oj] 


‘^Sxx 

”'^Sxy 

-Jc 

Sxz 

"'^Sxy 

'^Syy 

■'^Syz 

■'^Sxz 

"Jc 

Syz 

'^Szz 


(75c) 



C75d) 


where 


CjUsdZj 

"sy j 

Q ^s''s^^s 

r$z = 2 "ISL 5 

(76a) 

CjijdZj 

• 

“sy ' , 

1$ 

[(, »s"s<'7s 

• 

U/. =0 

Sz 

(76b) 


^ s 

“ '^Syy ' 3 "'s'-S ‘ L '’S*“s * ''$*‘*^5 

Ls Ls <760 

■'Sxy ° “ '^Sxz ° Jo ‘'syz ' jg ‘’s''s7s‘‘7s 
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f'-S 

^Sx " Jp Ps''s^s‘^^s 


h 


’Sy 




_ fS , 

^Sz " J« Ps^^s^'s " ''s‘^s^‘^^5 


(76d) 


in which p^ is the mass per unit length, is the length, and is the 
total mass of the uniform shaft. In addition, we can write the third term 

‘X. 

in Eq. (55) as 


{Ug} "" 2 


^0 




(77) 


Next, we wish to discretize Eqs. (76) and (77) and use the result to 
write a discretized shaft kinetic energy. Thus, let us assume that the 
shaft elastic displacements can be written as linear combinations of 
space dependent functions multiplied by time dependent generalized coordi- 
nates in the form 

''s ” *si(^s>M(t> <'®> 

♦s ' ®Si*^S*''*1**^* 

Note that «e are simulating the elastic shaft by degrees of freedom. 

We can write (nj) = figy ° ^"aSx *ASy “aSz^^ 

where w^Sx* '^ASy* '^ASz components of the vector w^^ along axes x^y^Zg 
in contrast to w^j^, w^g^, w^g^ components along axes x^y^z^. 

Substituting Eqs. (78) into Eqs. (76) and (77), we can write the shaft kinetic 
energy, Eq. (55), in the discrete form 
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" 2 '"S^'^aL '^ASy '^aL^ 6 '"S^S ■*■ “Sy^ 


s s 

I ^2 r r u’i’4’ 


*?"SZ I J, «Ij Kl%j * M\J> - “sx“sz j/r ■'ul 

■ %'"sz j, ^ I j, ”ij<^ui%j ^^1^o•> 

s s " 

^ '^ASx .l■^ '^ASy ^'^ASy^Sz “ ’^ASz^Sy^^J^ ^i’^ul 

^'^ASz^Sx ■ '^ASx^Sz^J^ ^Ki ^ I '"s‘-S^'^ASx“Sy ' '^ASy^Sx^ 


- n 


s s 

u u 


SX <"v1 ^ “Sy ,1, «r%< - "Sz J, jl, «n '^ul'-vj - %1%j> 

(79) 


where we introduced the definitions 

L. 


"ij - J 

M* = I 

' Jo 

L, 


0 "S^Si'^Sj'^^S 


M 


i^z _ 

i 


Ps^si^^^s 


'’s’‘'Si^S*^^S 


(80) 


Before we can write explicit expressions for the rotor blade "mass in- 
tegrals" Eqs. (61), we must apply the ordering scheme (46) to the blade 
kinetic energy Eq. (60). To this end, it is desirable to compare nondimen- 

sional quantities directly with the parameter e. Let us divide Eq. (60) by 
2 2 

the quantity mgfi Lg which has units of energy, where nig is the mass and Lg 
is the length of a rotor blade, and Q is the constant angular velocity im- 
parted to the shaft by the engine. Note that this division is compatible 
with the ordering scheme (46). 


34 
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We can now examine the rationale for retaining 0(c ) terms in the 
blade kinetic energy. This criterion is based on the derivation of Ref. 3 
which is concerned solely with rotor blades rotating uniformly about an axis 
fixed in inertial space. Hence, the derivation of Ref. 3 contains all of 
the third term and parts of the second and fifth terms in the blade kinetic 
energy, Eq. (60). Terms of order were retained while O(e^) terms were 
ignored. By examination of the fourth and fifth terms in Eq. (60), it is 
readily seen that O(e^) terms in {r^ + Ug) and {u^} can be of the same or- 
der of magnitude as the O(e^) terms of Ref. 3, depending on the magnitude 
of These terms are a direct result of the translational velocity of 

the blade root, i.e., of the point B, represented by the vector Wj^g in 
Eq. (43), and we note that in high speed flight can have large com- 
ponents, so that the fourth and fifth terms in Eq. (60) must include 0(e ) 
terms in {Ug} and {rg + Ug}, respectively. 

Substituting Eq. (48) into Eqs. (61), integrating over the undeformed 
blade, and retaining terms through O(e^), we can write the components of 
Eqs. (61) in a form analogous to Eqs. (76) for the shaft, as 


^By 

"B2 


p[x + u - e|^(v'+w’«>)cos - ejj^(w'-v'(>)sin ^gldxg 

p[6q + V + \ \ v'^)cos - e^(*+v'w')sin ♦gjdxg (81a) 

p[w + e^j^d - f } w'^)sin + e^* cos 
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'Bx 


^By- J 


'Bz 


p[u - e|^(v'+w'4.+w’*)cos (j.Q » e^(w'-v>-v'^)sin <j>Q]dXg 

p[v - e|^((j>^+v’v' )cos (j)Q - ej^(i+v'w'+v'w' )sin <l'Q]dXp (81b) 


p[w - e^( 4 -^+w'w')sin cos 4 >Q]dXg 


Bxx 


[p(®Q+v^+w^+2egv) + 2pe^(eQcos ♦q - e^^sin + vcos 


- v4»sin (>Q + wsin (^q + w^cos $q) + Jp]dXg 


^Byy 


f*"® 2 2 

^ [p(x +2XU+W ) + 2pe^(-xv'cos 4)g+xv'(jisin (fig - xw'sin ((.g 


xw 


2 

♦cos (|)g + wsin <fg + wKos ^q) + J 2 (sin ifig + 2({isin 4 >g cos ij)g 


2 2 2 2 2 
+ cos <J>Q - sin (|)g) + J^(cos 4>g - 2(f> sin (|ig cos 4>g 


+ (t^ sin^ 
L 


sin'" (|)g - cos^ 


'Bzz 


■B 


[p(x^+2xu+eg+v^+2egv) + 2pe^(-xv’cos *g + xv’^isin 4ig - xw'sin <fig 


- xw'4»cos <j)g + ®qCos <tig - eg<j»sin 4ig + vcos (j)g - V((.sin <t>g) 


0 0 

+ J 2 (cos'^ <^Q - 2^sin ^g cos 4>o - cos sin^ <^g) 


2 ...2 


. 2 


+ J^(sin" .(.g + 24 »sin .^g cos ^g - sin^ ^g + cos^ 
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'Bxy 


1 o 

[p(xeQ+xv+ueQ+uv) + pe^(xcos 4 »q - X(frsin (^q - ^ X(|. cos (f>Q 


1 uwl2. 


- 2 xv' cos <|)q - xv'w'sin <|)q - e^v'cos (p^ - eQw'sin (|>q 

2 

- vv'cos (()q - vw'sin <|>q + ucos <|>q) + pCj^(j)QSin (}.q + J2(-v'cos (|)q 

o 

- w'sin (j)Q cos (|)q) + J^(w'sin (|>q cos <J)q - v'sin <jiQ)]dXg 


(81c) 


Bxz 


•B 


12 12 

[p(xw+uw) + pe^(xsin (})q + x<j>cos (J>q - ^ X(j» sin <I'q - ^ xw' sin (})q 


- ww'sin (^q) - pC|^<l>oCos <Pq 


+ usin (j)Q - wv'cos (|>q 

J^(v'sin <pQ cos (Cq 


o 

+ J2(-v'sin (|)q cos <|)q - w'sin iJiq) + 


2 

w'cos <l>Q)]dXg 


'^Byz 


^B 

^ [p(eQW+vw) + pe^(eQSin (j)^ + eQ<l>cos ^ + vsin (^q + vKos (|)q 

0 


wcos <j>g - w<j)Sin (j>Q) + J2(sin (j)Q cos ((>q 

<j.g cos (j)Q + (|.s1n^ - *cos^ 


2 

+ <()C0 S «|»q) + J^(-sin <j>g 


((.Q - (frsin^ (|)g 


H, 


'Bx 


. . . . 

= J [p(wv-egw-vw)- Jp^ + pe|^(vsin + v^cos (^g - wcos (^g 
+ W(|)Sin (j)g - w^sin - Ogijcos tg - vKos <l>Q)]dXg 


H 


'By 


+ W(|)Sin (j)g 
^B 

[p(xw+uw-wu) 

0 


- w^sin (|.Q - 6g( 

xiji^sin (|)g - w'w'sin <j)g 


+ w 


w+uw-wu) + pe|j|{x|cos (|)g 

'wsin <|)g) + J2(v'sir 


+ ww'sin (Jig - usin (jig - v'wcos (jig - w'wsin (jig 

O 0 

+ w'sin'^ (jig) + J^{w'cos^ (^g - v'sin (jig cos <l>Q)]dXg 
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- f B . . . . 

^Bz " [p(eQij+vu-xv-uv) + p^(x^sin <f>Q + X(|i|cos <|>q - GqV'cos <|)q 

- eQw'sin <t>Q - vv'cos (|)q - vw'sin (|iq + ucos (|iq + xv'v'cos <|>j- 

Q 

+ xv'w'sin <1 )q + v'vcos (|)q + w'vsin ((iq) + J2(-v'cos‘^ ((iq 

• • 0 

- w'sin (|iq cos (|>q) + J^(w'sin (|)q cos (|)q - V'sin <I>q )]dxg 

in which we introduced the definitions 




2 

yc dnd^ 

^^2 = 

2 

yn dnd^ 

p = 


J 

Ja 



j, + = 


2 2 

y{n + r, )dndc pe„ 


A 

yndrid^ 


ydnd^ 


(82) 



yXCdndc 

A 


where y = y(x, n. C) is the mass density of the blade and the symbol A de- 
notes here integration over the blade cross-sectional area. Due to assump- 
tions of cross-sectional symmetry about the n axis and anti-symmetry of the 
warp function x, we also used the relations 


yCdndc = 

Ja 


yncdridc = 

'a 


JJ^ yXdrid? = 


Jj^ yXndndc = 0 


(83) 


In addition, we can write the third term in Eq. (60) in the form 

L, 


f ’B 

{Ub1^{Ub^^"’b ^ ^ 0 + 2pej^(-v^sin 


^ Wi^cos <('Q)]dXg 


(84) 


Equations (81) and (84) contain several O(e^) torsion terms which have 
been underlined. These terms are formally of higher-order and could be 
ignored. These terms were shown in Ref. 3, however, to be important for 
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low torsional natural frequencies, and hence, have been retained here for 
completeness. 

In actuality, to obtain Eqs. (81) and (84) we should integrate over 
the deformed blade or, equivalently, correct the integration over the unde- 


formed blade by introducing into the integrand the absolute value of the 

Jacobian associated with the transformation from the position of a point in 

the undeformed blade to the position of the same point in the deformed 

blade. The Jacobian in question differs from unity by a quantity of order 
2 

e , so that its introduction would not change our results substantially. 
Hence, the added complication appears unwarranted and will be left out. 

Next, we wish to discretize Eqs. (81) and use the result to writp t*" 
blade kinetic energy in a discrete form. To this end, let us conside» 
series expansions 




(85) 


where <ji ,, <}> . , . are space dependent admissible functions and q ., 

qyi» ^wi* '^(jii associated generalized coordinates. Note that we are re- 
presenting the elastic motion of each blade by N , + N + N + N^ degrees of 

U V W 9 

freedom. Substituting Eqs. (85) into Eqs. (81) and (84) we can write the 
discretized blade kinetic energy in a form resembling Eq. (79), where this 
form contains a very large number of integrals similar to the mass integrals 
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(80). The discretized blade kinetic energy is extremely lengthy and its 
presentation here is omitted for brevity. 

As mentioned earlier, a derivation of Lagrange's equations of motion 
requires also a derivation of the system potential energy. The potential 
energy is composed of two distinctly different parts; namely, the potential 
energy due to gravity, and the potential energy duetto the elastic effects. 
The following discussion is directed towards finding expressions for the 
potential energy of each component and ultimately the potential energy of 
the system. 

To calculate the gravitational potential energy of the airframe, we 
assume that the gravitational field is uniform. Hence, the potential energy 
of any mass point on the airframe is the product of the weight of the mass 
point multiplied by the height above the earth's surface, i.e., the verti- 
cal distance of the mass point from 0. This distance can be written in 
the form 

'’A = «0AZ * * “A> (86) 

SO that the gravitational potential is 


'^GA " ^ 




1 


= m^g(WoAz * 


[rA + u^ldm^) 


(87) 


where g is the local acceleration due to gravity. The last term of Eq. (87) 
is recognized as the vertical distance from A to the mass center of the de- 
formed airframe. When the elastic deformations are small, they can be 
ignored in Eq. (87) and one can replace the last term by the vertical dis- 
tance from A to the mass center of the undeformed airframe. 

The airframe elastic potential energy can be written in the general 


40 


form 


''eA ' 2 


( 88 ) 


where 

[K^] = [\][A^] (89) 

in which [a^] is the (diagonal matrix of natural frequencies squared, where 
the natural frequencies correspond to the elastic modes of the airframe. 

Of course, the matrix [K^] depends on whether free-free or cantilevered 
modes are used. Similarly, the matrix [M^] of generalized masses and the 
matrix [A^] of natural frequencies squared are different in each case. 

By analogy with Eq. (87), the gravitational potential energy of the 
shaft is 


Vgj = '“s9**'0AZ * ^^A0^^*'‘AS * “as*' * ^^SA^^ 




(rj , U5)dn,5 


(90) 


Assuming that the shaft undergoes normal stress along axis as well as 
shearing stress, so that letting EI^ be the flexural rigidity and GJ^ the 
torsional stiffness of the shaft, the elastic potential energy of the shaft 
is 



Substituting Eqs. (78) into Eq. (91), the elastic potential energy can be 


written in the discrete form 
V 


'ES 


-J I I kJ?" (n..,n 


]=1 jii iJJ VI vj 

in which we introduced the definitions 


S S 

+ + 2 A A K. . 

i=l j=l 


ij 




( 92 ) 
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(93a) 


, 0 ' 0 ' 

ij 



EIs^5^,".dZs 


(93b) 


where primes designate derivatives with respect to Zg. 

Assuming that the tail rotor is a rigid fan with its mass center at 
point T, the gravitational potential energy for the tail rotor can be written 
in the form 


VgT = "’t9(Woaz * ‘^AT^ ^ 

while the tail rotor elastic potential energy is zero by virtue of the rigid- 
ity assumption. Similarly, the gravitational potential energy of the rigid 
hub with mass center at H can be written in the form 




SH ■ 
(95) 


By analogy with Eqs. (87) and (90), the gravitational potential energy 
of a typical rotor blade can be written in the form 


^GB = "’b 9("0AZ " ^^3>"f^A0^'^»'AS " ^AS> " ^®3>[^A0^"t^SA^"^^SH " 


+ ^^FH^ ^^LF^^^®1^ 

B ^mg 

The elastic potential energy for a rotor blade is identical to that of 
Ref. 5 before the axial displacement u was eliminated. Assuming that the 
blade undergoes normal and shearing stresses the classical nonlinear strain- 

3 

displacement relations to 0(e ) are (see Ref. 3 ) 


42 


XX 


, ^ (v'2 + w'^) - x*” + * 2 ♦' ' 

- V[ncosC„ ^ ♦) - .sin(.o ^ *)] ‘ ^ ^ ^ 






For , linear stress-strain law. the blade elastic potential energy is si.ply 


'eB " 2 


Vcl 


[Ee^, t t c2^)]dnd;dx 


(98) 


«here E is the modulus of elasticity and G is the shear ™dolus. Substi- 
tuting Eqs. (97) into Eg. (98). integrating over the blade cross-section. 

and retaining terms as in Ref. 5 yields 


'EB " 2 


<E([U' t i (v'2 + W^)f A t W^.^■ + \ " ■>"% 

a v-^LI^cos^ ♦„ t I,sin2 . 2, cos *0(12 ' I,)] 

t w-^Cl^sin^ t I,cos2 ^2* sin ♦(, cos fodj ' I,)l 
+ 2[u' + f (v'^ + w’^)](*0*' + ^ 

. 2[u' t 1 (v'2 t w'2)][v"cos ♦„ * Wsin ♦„ - 2* sin 

- 2(+i+' t j 4.'^)[v"cos *0 t Wsin *,3 - 2* sin ♦(, cos ♦o(v"-W )]B2 
+ 2v"w"[sin (^Q cos <^q + <l>(cos <^q - sin ^qW 2 M 

- 2v"a“(sin ♦(, t Kos ♦o)C2 t 2WW(cos - fsin *0)02) f GIpW >dx 

(99) 


where A is the cross-sectional area and 


43 






2 . 2v2 


(n" + 
JA 


C ) dndi; 


'i = 


? dndr. Ip 

A ^ 


fl 


»2 = 


2 

Ti dnd^ 


n(n^ + C^)dndc 


Ae 


•f 

. = ndnd^ 

« JJA 


Ak 


■ II ..-II. I -gf 4 •£)']'* 


(100a) 


Cl = 


x^dndc 


= t^dndi 

JJa 


are certain cross-section integrals. Note that I-| and 1 2 are flapwise and 
chordwise area moments of inertia, respectively, Ak^ is the area polar 
moment of inertia. Ip is the torsional constant including cross-sectional 
warping and e^ is the tension offset from the elastic axis. We note that 
because of symmetry of the cross-section about the n axis and anti -symmetry 
of the warp function, we have taken into account in Eq. (99) the fact that 


• f ' 

Cdndr, = 

f 

ncdnd? = 

Ja 

Ja Jj 


;(n^ + C^)dndi; = 


rr 


xdnd; = 


Xndndc = 0 


Mnnhl 


Equation (99) can be discretized by using the expansions (85). Such a 
discrete expression for the blade elastic potential energy is very lengthy 
and will be omitted here for brevity. 

Equations (87), (88), (90), (91), (94), (95). (96), and (99) define 
the gravitational and elastic potential energy for all the system components 
In addition the shaft can undergo rigid body rotation relative to the air- 
frame, the flap link can rotate relative to the hub, and the lag link can 
rotate relative to the flap link. For generality, we shall assume that 
torsional springs are present to counteract these rotations. 

The potential energies due to these springs can be written as 


44 


'^KSx “ 2 '^Sk'^'x 

w K ,2 

'^KSy " 2 '^Sy^ 


'Kai 




'^Kei ' 2 

where k , k_ , k , k. are the torsional spring constants. 
Sx Sy a B 

system potential energy can be written in the form 


(lOla) 

(101b) 

(101c) 

(lOld) 
Thus, the total 


V ' ''ga * ''ea * ''gt * ''gs * ''es * ''gh * ''ksx ^ ''KSy 


* I " Vsi " 'gbi " ''ebi> 
1=1 


(102) 
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4. The Aerodynamic Loading and Virtual Work 


As indicated in Section 3, Lagrange's equations of motion require the 
expression for the nonconservative virtual work. In this section, we shall 
produce this expression for the main rotor, tail rotor, and airframe. 

We begin with the formulation of the main rotor aerodynamic loads. 

A strip theory or blade element approach is adopted because the results 
are in the form of loads per unit span, and thus are easily integrable over 
the total area of each blade, and because the theory allows a high level of 
flexibility in regard to the flow field assumptions. 

Before beginning the detailed analysis, a brief outline of assumptions 
and definitions is in order. The lift and pitching moment at a typical 
main rotor blade station are obtained via the Theodorsen expressions for os- 
cillating airfoils. The velocity used in these expressions is the local 
two-dimensional relative wind due to both nonoscillatory motions of the air- 
foil and induced flow, hereafter referred to as the nonoscillatory relative 
wind . In contrast, the blade element drag is assumed to depend upon and act 
along the total two-dimensional or oscillatory relative wind , which is defined 
as the velocity due to all blade motions and induced flow. In all phases of 
the aerodynamic formulations, effects due to reversed flow, stall, and shed 
wakes are ignored; compressibility is treated by use of the Prandtl-Glauert 
factor. 

In the induced flow analysis, it is convenient to work in terms of the 
hub x^^ plane rather than the average tip path plane so that the complica- 
tions of determining the latter are avoided. For the flight regimes to be 
studied, the orientations of the two planes should differ by only a few de- 
grees; this justification for using the hub plane remains to be verified. 
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The mean induced velocity, as determined from simple momentum theory 
(see Ref. 6) is 

W = Tp/(27.Lg^^B^') (103) 

where is the rotor thrust along the axis, Lg the rotor radius, the 
the air density, B the tip loss factor, and V the relative wind at the hub. 
The tip loss factor is given by 

B = 1 - -M- (104) 

where M is the number of rotor blades and 

Cf = 

The relative wind at the hub has the expression 

V = [(W-Wq^X ('^OAX 

where cij^ is the angle between the velocity of the hub and the plane normal 
to the hub axis z^. 

For hover, Wq^^^ = 0 and the induced velocity is 

= [Tj^/2TTY^B^Lg^]^^^ (107) 

For high forward speeds, W << Wq^^’ ~ '^OAX’ 

For moderate forward flight, Eqs. (103) and (106) are solved simultaneously 
for W, resulting in the quartic 
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( 109 ) 


^ sin * («0AX>^«^ - = 0 

The mean induced flow is then corrected by Glauert's expression 


“in ° • C “v *Z> U'O) 

o 

where r is the radial distance from the hub center and K , according to 
Payne (Ref. 7), has the form 


V - 


= r^ii.2 


x) 


( 111 ) 


in which 


M 


is the advance ratio and 


"OAX °R 


(112) 


X 


W, 


OAX 


sin 


“R 


fiL, 


(113) 


is the rotor inflow ratio. 

The induced velocity which is assumed to be in the -z^ direction, 
is then expressed in terms of components along the deformed blade section 
coordinate system x^y^z^ as 







(114) 


where [T^^^] is the matrix product [Tpg] [TgJ [T^^] (see Appendix A). 


The inertial velocity of a typical blade element is calculated by 
utilizing the previously introduced coordinate systems and velocity expres- 
sions. Specifically, the velocity of any point on the elastic axis of blade 
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i*, in terms of components along x^y^Zp, is given by Eq. (A33) of Appendix 

A with n and c set equal to zero. If this velocity is denoted by {wl)r. 

B E 

where the subscript E refers to the elastic axis, then the oscillatory 
relative wind (U) is given by 


{U} 





(W^)j 


('15) 


Since the profile drag is assumed to depend upon the two-dimensional 
oscillatory relative wind, we can use Eq. (115) to write its magnitude as 


^ ^ 2 ^DO^A^ ^ 

2 2 2 

where U = 1)2 + , c is the local blade chord length, and Cqq is the 

local profile drag coefficient (see Fig. 8). 

As mentioned at the beginning of this section, Theodorsen's equations 
call for the use of the nonoscillatory blade motion, which is defined as 
the sum of all motions except those which contribute to the unsteady velocity 
of the airfoil in the direction perpendicular to the nonoscillatory relative 
wind. This nonoscillatory motion is obtained by setting the flap angle rate 
B and the blade bending rate w (along with n and c) equal to zero in the 
right side of Eq. (A33) of Appendix A. If the resulting velocity is termed 
^'"*b'^EN ’ subscript EN refers to nonoscillatory motion of the 

elastic axis, then the nonoscillatory relative wind is given by 


{V} 



<«in' - 


'“b*en 


("7) 


As mentioned in Section 2, it is understood that the subscript i on variables 
associated with a particular rotor blade is deleted for clarity. 
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Theodorsen's expressions for the unsteady lift and pitching moment 


about the elastic axis are {see Ref. 8) 

L = C^^Y;^VbC(k)[Va^-i+b(l - a) a^] + 7rY^b2(Va^-Z-baa^) (118) 

and 

- [(^ - a) Vba^ + b^ (1 + a^) + baz] (119) 


In these expressions, is the local blade element lift coefficient per 

? 2 1/2 

unit length, V = (V2 ) ' , b is the blade semichord, a is the dis- 

tance from the midchord to the elastic axis in percent of semichord b, 
measured positive towards the trailing edge, C(k) is the complex lift defi- 
ciency function, is the instantaneous inclination of the chord to the non- 
oscillatory relative wind, and z is an upward (perpendicular to the nonoscil- 
latory relative wind) translational velocity of the elastic axis. The angle 
of attack can be written as (see Fig. 8) 


“A = 



( 120 ) 


The terms and z are given by 


“A = 


0 + 
cy 


cp 


( 121 ) 


z = 


J 

cos a. 


* ^DB31 


^ ^ ^DB32 '' ^DB33 


( 122 ) 


where B^ is the third component of the 3 x 1 matrix ([Tpp][rj-j^] [Tq^] x 

[ri^B^C^Lp] ’ '^bich is recognized as the coefficient 

of B in Eg. (A33) and T^g^i* ^0B32’ ^DB33 elements of the last row 

of matrix [Tpg]. The terms z and in Eqs. (118) and (119) will be 


neglected (see Ref. 8). 


KODUCmiUTY OF THE 
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The complex lift deficiency function C(k) is given by 


C(k) = F + iG 


(123) 


where 


and 


F 


(0,no>2 + (Y,-J„)2 


(124) 


G =- 


(J,+Yq)2 + (Y,-J(,)2 


(125) 


The Bessel functions , Jq, , and Yq depend on the reduced frequency 


k 



(126) 


where to, the frequency of oscillation is set equal to u for 1/rev oscilla- 
tions; more generally, to is the frequency which best represents the oscilla- 
tory motion of the main rotor blades. Because this frequency is not known 
in advance, an iterative process may be required, so that the value of (o 
assumed for the force computation is verified to agree with that resulting 
from the solution of the blade equations of motion. Note that C(k) intro- 
duces both a reduction in magnitude of the unsteady loads and a time lag of 


t 


lag 


tan 


6 

F 


(127) 


In the case of the quasi-steady assumption, C(k) = 1 and t^^^ = 0. 

At time t, the components of the aerodynamic forces in the y^Zp coor- 
dinate system are 
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( 128 ) 


(129) 


and the pitching moment is 

Re jMxoCt - 


(130) 


Now that the ai^rodynamic forces and moment acting on a typical blade 
section have been darived in terms of components along x^y^ZQ, we proceed 
to formulate the rnnconservative virtual work for the main rotor subsystem. 
The virtual work expression for the ith blade is 


6W = 


0 *''yD*"Dy ^d“oz * ”xD*°Dx> 


031) 


where 6Wpy and 6Wg^ are virtual displacements in the y^ and Zq directions 
and is a virtual angular displacement about the Xp axis. Explicit 
expressions for these quantities are listed in Appendix A as Eqs. (A39), 
(A40), and (A36). 

The fuselage aerodynamic characteristics are required for an eventual 
trim solution and for the airframe virtual work. The tail rotor exerts an 
aerodynamic force Tp on the airframe at point T. This force, which is 
assumed to act in the Zj direction (perpendicular to the tail rotor fan 
plane) and is calculated as part of the trim procedure, counte’^acts any 
airframe rotational tendency due to drag on the main rotor blades or other 
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aerodynamic effects. Also acting on the airframe at its center of mass 
are the overall airframe lift acting normal. to the distant free-stream 
velocity, the airframe drag acting along the distant free-stream velocity, 
and the airframe pitching moment about the y^ axis. These three airframe 
loads are given by 

“-A = 1 WAV%ef“AF (132) 

^A " 1 SA^AF^ef (133) 

^A " 2 ^MaAl^AW^^ref^SM^AF (134) 

where and are the airframe lift, drag and pitching moment 

coefficients based on the reference area is the static margin, 

a^p is the airframe angle of attack, and V^p2 = Due to 

the nature of the flight regimes to be studied, side forces are not con- 
sidered. 

The virtual work due to the three airframe loads and the tail rotor 
thrust is (see Fig. 9) 

‘“a ' <'-A **'■' “AF - "5 iftp) + 

(La cos sin o^p) ^ '*A«V 

V“aTz (135) 

where 6w^j^ is given by Eq. (A41) of Appendix A. 
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5. Lagrange's Equations of Motion in General Form 




I 

. 'JiP ’r.-' 


The kinetic energy, potential energy, and nonconservative virtual 
work have been discussed in detail in Secs. 2-4. In this section we shall 
present a general form of Lagrange's equations for the system. 

Before we can proceed with Lagrange's equations, it is necessary to 
identify the functional dependence of T and V on the various generalized 


coordinates. The functional dependence presented here is the most general 

possible and not restricted to the linear case. First, we recognize that 

the airframe is described by three rigid-body translations w , w , w 

OAX OAY OAZ 

three rigid-body rotations x^, x^, and P elastic generalized coordinates 
(i = 1,2,. ..,P). Then, considering Eqs. {65a, b), (69), (70), (71), (87), 
and (88) in conjunction with Eqs. (3) and (4) as well as the definition of 
we conclude that the functional dependence of T., Vp., and Vr-. is as 

n bn t A 

follows: 


^A ^A^'^0AX»'^0AY '^OAZ’^x’V’^z’'^Ai’^x’^y’^z^ ’ i = 1 .2 P 


''gA ' '^GA^'^OAZ’^x’V’^z’^Ai^ ’ i ~ 1’2,...,P 
, i = 1,2,.. .,P 


^EA = 


(136a) 

(136b) 

(136c) 


The transmission shaft is described by three rigid-body rotations 
ip , and ip relative to the deformed airframe and 2S + S elastic generalized 
coordinates (j = 1 ,2, . . . ,5^^; k = 1,2,.. .,S^). Note that the 

rotation = Qdt is a specified function of time. Considering Eqs. (79), 

(90), (92), and (101a, b) for the shaft kinetic energy, gravitational poten- 
tial energy, elastic potential energy, and potential energies due to torsional 
springs at the base, respectively, together with Eqs. (3), (4), (11), (12), 
(20), (69), (78), (A14), and (A30), we can write the functional dependence 
^S’ ViS’ '^ES’ '^KSx '^KSy follows: 
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's " 's^'^OAX*'^OAY* OAZ’ x’V’^z*'’Ai’'*'x’V'’uj’%j’ X* 

i = 1,2....,P; j = 1,2,...,S„ (137a) 


''gS = '^GS%AZ»\*V’^z’%-^’^’\j’"vj^ 


1 ~ l»2j...jPj j — l>2(...tS, 


(137b) 


% ^ " 1,2,...,$^; k = 1,2,...,$^ 


'^KSx ^ \sx^*^x^ 


(137c) 


(137d) 


V = V., (il» ) (137e) 

TSy Ksy^%^ 

The tail rotor is described by one specified rotation | dt relative 
to the deformed airframe. The kinetic energy of the tail rotor is given by 
Eq. (57), so that considering Eqs. (3), (4), (26), (27), (28), and (66) the 
functional dependence of Tj can be written as 

I 

"•^T " ^T^'^OAX’'^OAY*'^OAZ’^x’^y’^z’’^Ai’^x’^y*^z’’^Ai^’ " ‘•>2,...,) (138) 

On the other hand, the tail rotor potential energy is given by Eq. (94), and 


it is readily seen that 

'^GT ^ '^GT^'^OAZ’^x’^y’^z’^Ai^ ’ ^ " ^*2,...,P 


(139) 


Following the same pattern, the hub is fixed to the end of the transmission 
shaft and moves together with it, so that considering Eqs. (58) and (95) in 
conjunction with Eqs. (3), (4), (11), (12), (20), (21), (69), (78), (A15), 
and (A31), we conclude that the functional dependence of T^ and Vg^ is 

^ ^H^'^OAX’^OAY^OAZ’^x’^y’^z’^Ai’'^x’'^y’'^uj’'^vj’%K’\’V’^z’^Ai* 

'l'x’’^y »'^uj ’*^vj *'^^k ^ ^ ” 1»2,...,P; j = 1,2,. ..,S^; 


k = 1 ,2, . . . ,S. 


(140a) 
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'^GH " ^ " 1.2,...,P; j = 1,2 

(140b) 


Finally, the rth rotor blade is described by the three rotations 
(lead-lag), (flapping), and 9^ (pitch) as well as the % 
elastic generalized coordinates q^^^, q^^^ {a = l,2,...,N^j; 

m = 1,2,...,N ; n = 1,2,...,N s = 1,2,...,N.). As mentioned in Sec. 2, 

V W 9 

of the three rotations only and are generalized coordinates. The 
kinetic energy is given in general form by Eq. (60). Considering Eqs. (8i), 
(84), and (85) along with Eqs. (3), (4), (11), (12), (20), (21), (30), (31), 
(36), (40), (69), (78), (A16), and (A32), the functional dependence of 
Tgr (>" = 1,2,...,M) is deduced to be 


^Br " '^Br^'^OAX’'^OAY’Vz’\’V’^z’'^Ai’'‘^x’'^y’%j’\j’V’“r’®r’V£’%rm’ 

• • 

*^wrn’'^(|irs’^x’^y’^z’'^Ai *'^x’'^y’’^z’'^uj ’"^vj ’'^4)k*“r’®r’^ur^’^vrm’ 
\rn’‘^(^rs^ ’ ^ ^ 1*2,...,P; j = 1 ,2, . . . ,S^^; k = 1,2,...,$^; 


S, - 1j29t*«iN » ni “ 1y2»«.»9N 9 n - 1y29.*.«N « S - 1)2y».»yN, 

U V W (p 


(141) 


The gravitational and elastic potential energies Vgg^ and are given by 
Eqs. (96) and (99) respectively, whereas the potential energies and 
'^Kpr given by Eqs. (101c) and (lOld). Their functional dependence is 
readily seen to be 


“ (w_. 7,X„,X„,X,,n/ii »n,,i,n..j,n..,a^,P„,q 


GBr ’'GBr^'^OAZ’ x’ V z’''Ai’^x’'^y’''uj’"vj’"^k’“r’'’r*\r)l’'^vrm’\rn’>rs 


» 


i = 1,2 P; j = 1,2 S, ; k = 1,2,. ..,S^; I = 1,2, ...,N„; 

u $ u 

m - l,2,...,N^i n - 1,2,... »Ny^» s - 1,2,. ..,N^ 


(142a) 
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y "Vf«(Q *Q >Q >Qi )iy-'“l*2*.**jN>m“l>2f 
^EBr EBr^^uri ^vrm* wn*^(|»rs' u 

n = 1»2,...»N^» s “ 1»2»...»N 


•>N,; 


<t> 


'^Kar " '^Kar^“r^ 


'^Kur " '^K3r^^r^ 


(142b) 

(142c) 

(142d) 


Considering Eqs. (62) and (102), the system Lagrangian is 

M 

L = T-V = h^h^h^V\^-h^-\r\s-hs-\H-\sy.-\sy ^ 
■'^6Br"'^EBr''^Kcxr"'^KBr^ 

Noting the functional dependencies (136) - (142), Lagrange's equations for 


the 

three 

rigid-body translations can be written as 
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3L 

d 3 

<Ta * - Tj . T„ f 

^Br> = 

dt 

3Wqay 

SWqay 
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3 Wqay 
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3L 

3L 

_ d 3 

<Ta *h*h^^H* j, 


dt 

^WqaZ 

^'^OAZ 

dt „• 

^wqaz 


(^ 3 , . V 3 , . V 33 . V 3 , . V 33 ^) = 


(144c) 


where Fq^j^, Fg^y, Fq^^ nonconservative generalized forces associated 

with Wq^j^, WQ^y, Wq ^2 f'espectively. Similarly, we can write the rotation 
equation 

M 


1_ ik 

n. 


3l . d 3 

n 


^ . Tj ^ T„ . T3,) - 


[T, + T, 


3X^-A 'T 


M 


+ T$ + Th - VgA - - Vgj - Vg(^ * (Tj^ - Vgg,,)] - (145) 


57 


where analogous equations can be written for \ and X and F , F. , 

y z XX Ay 

are the nonconservative generalized forces associated with X , X , X 

X y 

respectively. In the same manner, we can write the airframe elastic 
mation equations 


dt 


3L 



8L _ 3 

3HAi ■ dt 3nAT 



” 9HAi ^ ■ ''gA " V ■ ''gS ■ '^GH 

' '^EA ^^Br " '^GBr^^ ^ *^nAi ’ ^ 

r= I 

the transmission shaft rotation equation 


d 3L 3L _ d 3 
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"t ■ « si * X ^Br> - 9*^ - \s 

X 

M 


'KSx ^^Br ■ '^GBr^^ ^ 


where an analogous equation can be written for the transmission 
elastic deformation equations 
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z’ 
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(147) 

shaft 

(148a) 

(148b) 
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1 _ 9 L 
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%k - « ■ r ^1 "»'• 


■ ’'eS^ ' ''nH 


j ” l»2»...»Syj k - 1t2»...iS 


the rotor blade rotation equations 
3L 
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and the rotor blade elastic deformation equations 
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r = 1.2... .,M 

where the meahing of F^^. F^^j, F^^j. F^^^, F^^. F^^, 

^qvrm’ ^qwrn’ '^qtrs “ 1.2.. ...P; j = 1.2 S^i k = 1.2 S^; 

t ~ l»2»...*Nj m = 1»2 j..4,Nj n = l,2t...tNi s - l|2t...,Ni r - 
u V w (> 

1,2,... ,M) is obvious. 


Equations (144) - (150) represent a set of 6 + P + 2 + 2S + + (N 

U u 


59 


+ Ny + Nw + N^)M coupled nonlinear equations which must be solved simul- 
taneously. Before we proceed with the actual solution, we must first ren- 
der the equations in a more explicit form by replacing the various kinetic 
energy components, potential energy components, and generalized noncon- 
servative forces by their specific matrix expressions. Then, we derive the 
variational equations by expanding the nonlinear equations about trim 
solutions. A general form of the trim solution and variational equations 
will be discussed in the next section. 

Lagrange's equations, Eqs. (144) - (150) Involve an extremely large 
number of matrix multiplications and differentiations both with respect 
to generalized coordinates and velocities and time. The magnitude of the 
task demands a more automated approach, so that in Sec. 7 we present a 
procedure for the derivation of Lagrange's equations by computer manipu- 
lation. 
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6. The Perturbation Equations 


Let us assume that, following discretization, the helicopter can be 
represented by an n-degree-of- freedom system, so that its motion is des- 
cribed by n second-order Lagrange's equations or 2n first-order Hamilton's 
equations. The latter set can be written in the form 

's. 

^i " Xi(><l*X 2 ,...,X 2 ^,t) , i = 1,2 2n (151) 

where are generally nonlinear functions of the variables x. (i = 1,2,..., 
2n) and of the time t. Note that n of the variables x^. represent generalized 
displacements and the remaining n represent generalized velocities, or 
generalized momenta. The 2n variables x.(t) define the state of the system 
at any time t. 

Next, let us consider a special solution of Eqs. (151), namely, a 
trim solution. In general, there are many such solutions, but solutions of 
particular interest are those corresponding to hovering and to forward 
flight. Denoting a particular trim solution by (t), and recognizing that 
<|>^(t) must satisfy Eqs. (151), we can write 

<l>^(t) - X. ((f)^ » i “ 1,2,. ..,2n (152) 

In general, such solutions are periodic, iti^.(t) = (j)^(t + T). We shall refer 
to as the unperturbed motion . 

Letting y^(t) be perturbations about a given trim solution <))^(t), the 
general perturbed motion can be written in the form 

x^(t) = 4>.(t) + y.(t) , i = l,2,...,2n (153) 

so that, introducing Eqs. (153) into Eqs. (151), we obtain 
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^ y}*^2 * ^2’*‘*’^2n '*’ ^2t. ^ ' l»2....,2n (154) 

Considering Eqs. (152), Eqs. (154) reduced to 

“ X^.((^i + yi»^2 ^ ^2’‘‘*’^2n ^ ^2n’^^ ” X^ (^^ ,4>2» • • • »42p*^^ 

i = l,2,...,2n (155) 

which are referred to as the differential equations of the perturbed motion . 
Equations (155) can be expressed in a different form. To this end, let us 
expand the first term on the right side of Eqs. (155) in the Taylor's series 
about the trim solution 


X^ ('t'l yp'i*2 ^ ^2’***’*^2n ^ ^2n’^^ ~ X^ ( 41 ^ ,^2» • • • 


2n 3X 


^ .1 W 
j=l ""j 


x=$ y^- + Oi(r) . i = l,2,...,2n 


(156) 


where x, and y are 2n-dimensional vectors associated with x^, 4>^, and y^ 

2 

respectively, and 0^(y ) denotes terms of second order in y^ . Introducing 
the notation 


3X 




, i,j “ 1,2,..., 2n 


(157) 


and considering Eqs. (156), Eqs. (155) can be rewritten in the form 
2n , 

y^ = I aij(t)yj + ) • 1 = 1.2....,2n 

J " 1 


(158) 


where in general the coefficients a^j are periodic, a.j(t) - a^j(t + T). 

Note that Eqs. (158) are nonlinear because of the terms 0.(y ). 

A case of particular interest is that in which the perturbations y^(t) 
are small. In this case, we can neglect the second-order terms in Eqs. (158) 
and obtain the set of linearized equations 
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(159) 


2n 

^ ~ l»2,...»2n 

which are referred to as the variational equations . 

The perturbation equations, Eqs. (158), or the linearized version, Eqs. 
(159), were derived on the assumption that the solution (|)^. (t) represents an 
actual trim solution, i.e., they represent a solution of the original equa- 
tions (151). Trim solutions, however, are difficult to obtain and at times 
one may wish to assume an approximate solution and derive a set of pertur- 
bation equations about the "assumed trim". The question arises naturally as 
to the behavior of these equations. To answer this question, let us denote 
the actual trim by <j)^(t) and the assumed trim by <}i^(t). Assuming that the 
two solutions differ to some extent, we can write 


♦J(t) = *,( 1 ) ♦ «,{t) , 1=1.2 2n (160) 

where «.(t) (1 = 1.2 2n) represents the difference between the two sole- 

tions. Then, the perturbed motion can be written in the form 


x.(t) = t*(t) + y*(t) , i = l,2,...,2n (161) 

where y*(t) are perturbations from the assumed trim. Inserting Eqs. (161) 
into Eqs. (151), we obtain 

♦*(t) + y*(t) = X^(4>* ^ y ^,<>2 + y2’--***2n * ^2n’^^ 

i = 1,2,. ..,2n (162) 

so that, expanding about the assumed trim we can write Eqs. (162) in 
the form 

2n 3X 

y^-(t) = X^(it)^,(Ji 2 f Ix^ x=^* ^i * ®i^- ^ ’ 

i = l,2,...,2n (163) 
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Unlike the case in which the expansion was about the actual trim, however, 
the first two terms on the right side of Eqs. (165) do not cancel out, be- 
cause the assumed trim does not solve Eqs. (151). It will prove of interest 
to examine these terms a little closer. 

Assuming that the difference between the actual and the assumed trim is 
relatively small, we can write the expansion about the actual trim 
* * * 2n 

(^l ,^2 ’ • • * ~ ’*^2’ * * *^2n’^^ ^ 

i = 1,2,. ..,2n (154) 

where 


3X, 


' 3x, 


x=4> 


(165) 


are the actual coefficients, which are g eneral ly not known . Considering 
Eqs. (161) and (164), as well as Eqs. (152), we obtain 


* * * .* 2n 

X^‘ (^l t4>2» • * • “ ^^(^) ~ I ^ik^^^^k^^^ ” 

k“ 1 

Moreover, introducing the notation 
2n 

£^(t) " ^ ” 4^(t) , i - l,2,...,2n 


(166) 


* 9X. 


(167) 


Eqs. (163) reduce to 

y*(t) = I 
3 1 

Hence, e^(t) play the role of unknown extraneous forces introduced by the 
process of using an approximate trim Instead of an actual one. Although 


(168) 
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these forces are not generally known, one may be able to estimate them. 

•k 

The solution y^. represents perturbations from the assumed trim instead 

of from an actual trim. The question can be asked as to how they compare 

with the perturbations y. from the actual trim. For small deviations 6. 

★ 

from the actual trim, the response y^ should not differ very much from y^. , 

★ 

but this cannot be taken for granted. The relation between y^. and y. de- 

* 

pends of course on Methods for estimating bounds for y - y for given 

•k 

bounds for $ - appear highly desirable. 
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7. Algebraic Computer Manipulation 


In Sec. 3 we mentioned two approaches for expanding the system kinetic 
energy, Eq. (63). One approach is to obtain explicit expressions for the 
mass integrals given by Eqs. (53), (56), and (61) while retaining the trans- 
lational velocity and angular velocity of each set-^of axes in implicit form. 
Such an approach is convenient when the elastic displacements are of inter- 
est and was the approach used in Sec. 3. The other possibility is to sub- 
stitute explicit expressions for the translational velocity and angular 
velocity of each set of axes given by Eqs. (A28)-(A32) and (A12)-(A17) into 
‘•qs. (64) while retaining the mass integrals (53), (56), and (61) in im- 
plicit form. This approach is convenient when the rotational coordinates 
as well as the terms due to coupling between bodies are of interest, and it 
must be adopted if one wishes to derive explicit expressions for Lagrange's 
equations. Examining Eqs. (A28)-(A32) and (A12)-(A17) one concludes inme- 
diately that working with explicit expressions for the velocities and angu- 
lar velocities of each set of axes involves the calculation of very lengthy 
matrix products. Moreover, these matrix products involve quantities of 
different orders of magnitude and/or importance. Many of these terms are 
insignificantly small and can be ignored. Because of the complexity a,sso- 
ciated with algebraic multiplication of a large number of matrices and be- 
cause of the high probability of human error in performing these multipli- 
cations by hand, it seems highly desirable to computerize these algebraic 
operations. It should be pointed out that algebraic computer manipulation 
need not be restricted to matrix multiplication. Indeed, many operations 
involved in the derivation of Lagrange's equations, including differentia- 
tion, can be performed by computer. 
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nEPRODUCmmiTY OP THE 
ORIGINAL PAGE IB POOR 
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In this section, we shall outline an algebraic computer procedure which 
can be used to expand the system kinetic energy explicitly in terms of the 
translational and angular velocity of each set of axes in such a way that 
small terms are ignored automatically. In addition, the same method can be 
applied to the gravitational potential energies, Eqs. (81), (90), (94), 

(95), and (96). Furthermore, the system kinetic energy and gravitational 
potential energy expressions obtained by the computer method are in a form 
which can be easily differentiated algebraically by the computer, so that 
ultimately Lagrange's equations, Eqs. (144)-(150), can be obtained explic- 
itly with a minimum of human effort. 

The ideas used in implementing algebraic manipulation on a computer, 
are best introduced via an explicit example. Let us consider the product 

-0.5 6 se Ccx(-6 Sa - 6 C6 Ca - B Jg^^ 

which appears in the expansion of {Og}^[Jg]{fig) , where the latter was 
encountered in Sec. 3. To calculate this product on the computer, we shall 
associate numbers with the different groups of characters that represent 
quantities to be manipulated algebraically, i.e., B, ‘^Byz’ '^Bzz’ 

c0, sa, ar ca. These groups of characters constitute what will be re- 
ferred to as symbols . The association of numbers with symbols allows us 

I 

to substitute manipulation of numbers for manipulation of symbols and is 
accomplished by the formation of a symbol table . Table I. The symbol table 
contains two entries per line. These entries give the character representa- 
tion of a symbol and a weight assigned to the symbol, where the weight is 
determined by the analyst according to his knowledge of the symbol's mag- 
nitude or his desire to retain its effects. The weight need not be a fixed 
quantity and can be changed at will. A high numerical value of the weight 
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implies a high-order term, i.e., a loss significant term. Note that here 
we assigned so, co, and ca a weight of zero, sa, and 8 a weight 

of unity, and a weight of two, where the weights represent anticipated 
magnitudes of these symbols. The number associated with a particular sym- 
bol is the line number in the symbol table. In this case, one is associated 


with 6, two with with Jg^^’ 

Examining Eq. (169), we see that it is necessary to form algebraic 
products of symbols, for example - 0.5 8 se ca. To this end, we define 
terms . A term consists of a signed numerical coefficient , a pattern con- 
sisting of the numbers associated with each symbol appearing in the product, 
and a weight which is the sum of the weights of each individual symbol ap- 
pearing in the term. Hence, in view of the number-symbol associations and 
the weights of Table I, we represent -0.5 6 S9 ca as a term having a coeffi- 
cient of -0.5, a pattern of 1 , 5, 8, and a weight of 1 . All terms, i.e., 
all coefficients, weights, and patterns are stored in numbered storage 
stacks. The coefficient and weight of any term are always single numbers 
and are stored side-by-side in the coefficients and weights free storage 
stack , exhibited in the form of Table II. On the other hand, the pattern of 
a term may differ from the pattern of other terms and must be able to re- 
present the product of any number of symbols. Because of the different 
lengths of different patterns, all patterns are stored in the separate 
patterns free storage stack , labeled as Table III. Note that the coefficient 
and weight of -0.5 8 se ca are stored in line 3 of Table II and the pattern 
is stored in lines 5, 6, and 7 of Table III. 

It is also necessary to form algebraic sums of terms, such as -8 Jg^^ sa 
- 6 Jgy 2 c0 Ca - 6 Jg ^2 se Ca, which we shall call series . Each series is 
given a distinct series name . As examples, we shall call the series con- 
sisting of the single term -0.5 8 se ca by the name X, the series -e Jd sa 

dXZ 
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- B J„ CO Cci - e Jd s0 ca by tho name Y and the series resulting from 

Byz bzz 

the product of X and Y by the name Z. A series is described by a sequence 
of terms with coefficients and weights stored sequentially in the coeffi- 
cients and weights free storage stack and with patterns stored sequentially 
in the patterns free storage stack. The summation of terms in the sequence 
is understood. Hence, the series Y is described by the coefficients and 
weights stored in lines 8, 9, and 10 of Table II and the patterns stored in 
lines 12-14, 16-19, and 21-24 of Table III. To distinguish sequences of terms 
forming series, each series name is assigned a number corresponding to a 
line in the series definition table , Table IV. Each line of this table 
contains three entries giving the line number of the coefficient and weight 
in the coefficients and weights free storage stack of the first term in the 
series, the line number of the beginning of the pattern in the patterns 
free storage stack of the first term in the series, and the number of terms 
in the series. Assigning the number 5 to the series name Y, the fifth line 
of Table IV contains the entries 8 and 12 giving the storage locations of 
the series and the entry 3 designating that there are three terms in the 
series . 

Let us now consider the multiplication of two series, namely, the mul- 
tiplication of Y by X, which can be performed term by term. The product 
of two terms yields a new term. If the total weight of the new term, given 
by adding up the weights of the two terms in the product, is greater than a 
specified value, for example 3, then the new term is deleted. Otherwise, the 
coefficient of the new term is the product of coefficients and the pattern 
of the new term is the concatenation of the patterns of the two terms in the 
product. Multiplying tl’e first two terms in Y by the single term in X, the 
total weight of each resulting new term is 4 which is greater than 3, so 
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that these terms are deleted. Multiplying the last term in Y by the single 
term in X, the resulting new term has a total weight of 2, a coefficient of 
0.5, and a pattern of 1, 1, 4, 5, 5, 8, 8. Assigning the number 6 to the 
series name Z, this new term which is the product of X and Y is stored 
according to the information in line 6 of Table IV. 

The method outlined can be programmed easily in Fortran IV and is 
appealing because of its simplicity. We have discussed only mulitpl i cation 
of series. Clearly, considering each entry in a matrix to be a series, al- 
gebraic multiplication of matrices is accomplished by multiplying and adding 
series. In addition it is not hard to see that differentiation is simply 
a matter of looking for the occurrence of particular symbols in the pattern 
of each term. In future work, we shall present detailed documentation of 
a computer program which performs algebraic manipulation and we shall use 
the program to obtain explicit expressions for Lagrange's equations. 
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Table III - Pattern Free Storage Stack 
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Y = 5 
Z = 6 
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Table IV - Series Definition Table 
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8. Summary and Future Plans 


This report presents a formulation of the equations of motion of a heli- 
copter. The method of approach is a variation of the component-mode synthesis 
in the sense that it regards the aircraft as an assemblage of interconnected 
substructures. The substructures identified are the airframe, the trans- 
mission shaft, the tail rotor, and the main rotor. The rotor blades are 
assumed to be articulated with the flap-lag-pitch configuration. The equa- 
tions of motion are derived in general form by means of the Lagrangian 
formulation in conjunction with an orderly kinematical procedure that takes 
into account the superposition of motion of various substructures, thus 
circumventing constraint problems. 

Because of the complexity of the problem, the derivation of explicit 
equations of motion is sure to be extremely tedious and time consuming. 
Moreover, the probability of error in deriving the equations is large indeed. 
Fortunately, a number of assumptions can be made to simplify these equations. 
In particular, one can ignore certain higher-order quantities in the equa- 
tions for the rotor blades. However, this task is also sure to be tedious 
and time consuming. Hence, a procedure for the automation of the derivation 
of the equations of motion is unavoidable if time and effort are to be mini- 
mized. Such a procedure consists of a computer program capable of performing 
the many matrix multiplications involved, certain differentiations, and elimi- 
nation of higher-order terms. This latter task can be made easier by 
adopting an ordering scheme. In a computer manipulation the ordering scheme 
can be altered, thus producing sets of equations corresponding to different 
sets of assumptions. 

The next phase of the investigation is concerned with the derivation of 
the equations of motion in explicit form. To this end, a method for the 
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derivation of the perturbation equations (which are generally nonlinear) by 
means of computer manipulation will be developed along the lines of Sec. 7. 

Implicit in the derivation of the explicit equations of motion is 
the truncation problem. Truncation will be done first on the substructure 
level and then on the complete aircraft level. An important question is 
that of the airframe modes and how they can be used to evaluate certain "mass 
integrals." In this regard, it may prove advantageous to look into the 
possibility of using "admissible vectors" to represent the motion of the air- 
frame instead of using actual airframe modes. 

Another problem that needs to be answered is the effect of nonlinear 
terms in the perturbation equations for the blade motion. If such terms 
cannot be ignored, considerable difficulty is likely to be encountered in the 
determination of the dynamic characteristics of the aircraft. Intimately 
related is the question of the axial displacement of the helicopter blade, 
as the axial displacement introduces nonlinear terms. The question is 
whether one should treat this axial displacement as an independent distrib- 
uted coordinate or attempt to express it in terms of the bending displace- 
ments. Of course, if the axial displacement can be expressed in terms of the 
bending displacements, then this fact in itself implies a reduction in the 
number of degrees of freedom of the simulation, as the differential equation 
for the axial displacement is eliminated for every blade. Under certain 
circumstances, it may be possible to eliminate also the differential equation 
for the torsional motion in a similar fashion. 

It appears desirable, if at all possible, to derive a set of linear 
ordinary differential equations with constant coefficients for the system, 
as such a set leads to an eigenvalue problem likely to yield useful informa- 
tion concerning the helicopter dynamic characteristics, such as natural 
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frequencies and natural modes of vibration. Then, it may be possible to 
treat nonlinear effects as perturbations on the linear case. If a set of 
linear differential equations with constant coefficients can be obtained 
for the system, then it must by necessity correspond to the case of hovering. 
In the absence of aerodynamic forces, the set of equations is bound to be of 
gyroscopic type. The solution of the eigenvalue problem for the gyroscopic 
system can be obtained by the method of Ref. 9. The possibility of using 
the natural modes of the gyroscopic system to truncate the overall problem 
will be explored. 
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re 1. The Fuselage Generalized Coordinates 



Figure 2. The Fuselage Coordinate System 
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^ COlRfeCTtftN 6F TftlU 


Rotor, thrust) 

Figure 5. The Tail Rotor Coordinate System 
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Figure 7. The Deformed Main Rotor Blade 
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Figure 6a. The Main Rotor Hub and Flap 
Coordinate Systems 


6b. The Lag Coordinate 
System 


6c. The Undeformed Blade 
Coordinate System 




Contained herein is a summary of the key translational and angular 
velocity relationships developed in Sec. 2. Also included, when needed 




elsewhere in the text, are the full expansions of these expressions in terms 
of generalized coordinates, control variables, and elastic displacements. 
Important virtual displacement expressions are also listed. As mentioned 
in Sec. 2, it is understood that the subscript i on variables associated with 
a particular main rotor blade is deleted for clarity. Because many of the 
matrix products appear repeatedly, the following compact notation is de- 


veloped; 

[Top] = 

The important angular velocity expressions in Sec. 2 are: 
i^f^} = X^EXylCx^Ke^} + x^[Xy]{e^} + (A3) 

+ [%][v]([ilQj]{u^j}) + 

-v^ (Lj,t)^ 

{51h> = ^ E4»t) / (A6) 

, ELg.t)) 

{fip} = [Tp^lt^H^ ' 

i\} = [T^pli^lp) + a{ 63 } (A8) 
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o(e^ 

} 







(A9) 


+ 

(u^d) 








(A10) 

w'v' 

+ ^ 










(v- - 

;v^ ■ 


'w' 

2 

)sin (|»Q - 

(w- 

- vV 

'w' 

- ;■ 

- ]■ w'v 

2 

'* )cos ((.q 

.(v‘ - 

vV^ • 

■ io 

•w* 

^)cos .Jq + 

(w- 

- vV 

'w' 

- w' 

,2 1 •, 

W - ^ W V 

'■^)sin (|.g 

+ 

v' ((j)COS 

^0 - 

1 

2 

sin (|)q) 

+ w' 

(({isin 

^0 

4 

2 

<|> cos (|)q) 

' (All) 

- 

v' ((jisin 

^0 " 

1 

2 

cos (j)^) 

+ w' 

' (i|>COS 

^0 

1 

■ 2 

2 

<t> sin ^q)^ 



Substitution of Eq. (A3) into Eq. (A4) yields a fully expanded expression 
for {n-p}, substitution of Eq. (A3) into Eq. (A5) yields a fully expanded 
expression for etc. Repeating Eq. (A3), this process yields the 

following equations: 

(f^A^ = xJXy][x^]{e3} + + Xy{e2> (A12) 

{Qj} = X^[Tp^][Xy][x^]{e3> + X^[T3-^][x^]{e^} + ^y[T-p^]te2) 

+ [TjA][v]([tgy]{u^x}) + 

{fig} = 

+ [T 5 l^][v]([^Q 5 ]{u^ 5 }) + } + iJiyEi^^Ke^} + ^{ 03 } (A14) 

{n^} = X^[T^^][Xy][x^]{e3) * ■*■ \*-'''hA^^®2^ 

ps 

+ il[T|^s](e3} + [Tpij] I U5 ^*-5’^^ / 
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+ [Tg^][v]([£gs]fjAS^^ ^ ^ 

+ ^[Tg5]{e2) + [Tgg] < uj (Lg,t)|- 6[Tgp]{e2} + “[TgLlte^} 

\ h X 

+ e{e^} I 

-vj (Lg.t) 

+ ^[Tpglfe^} + [Tpg] , Ug (L^,t) ' - B[Tpp]{e2} + a[TgL^^®3^ 

1$ (Lg.t) 


(A16) 


+ 6[Tgg]{e^} + yv' - v'v'^ - ^ v'w'^)s1n 

((V' - VV^ - |- V'W'^)C 0 S (t»Q 

w'v' + 4. 

- (w* - vVw' - w'w"^ - ^w'v"^)cos 4p + v'Ckos 4q 

- (w* - v'vV - w'w'^ - ^wV^)sin 4^ - v'(4sin 4^ 


12 1 ? 

- 2 4 sin 4 q) + w'( 4 sin ^q + ^ 4 ^ cos 4q) 

+ I 4^ cos 4q) + w'( 4 cos 4q - ^ 4 ^ sin 4 q) 


The important translational velocity expressions in Sec. 2 are: 

{WoA> = [Wqax Wq^y '^OAZ^^ 

{Wa> = ‘ ^''a ^ ‘^A^^‘"‘A^ * ^^A^ 


(A17) 


(A18) 


(A19) 


ROTODUCmiLlTY OF^ 

original page IB POCHt 




^ ^^AO^^'^OA^ ■ ^*"aS ^‘^AS 


{WsH> ^''’sA^^'^AS^ ■ ^*"SH * ^ ^^SH 

{Whp> = t'^'HS^^'^SH^ ■ 


{Wp^} = [Tp^]{w^P> - [rp|_]{np} 




{Wg} = [Tg|^]{W|^g} - [fg + Ug]{Qg} + {Ug} 


} 

} ■• 
} 


(A20) 

(A21) 

(A22) 

(A23) 

(A24) 

(A25) 

(A26) 

(A27) 


Substitution of (A18) into (A21) yields a fully expanded expression for 
{w^g}, substitution of the results into Eq. (A23) yields a fully expanded 
expression for iwg^}, etc. Note that the use of Eqs. (A3)-(A9) is necessary. 
Repeating Eq. (18), this process gives the following selected equations: 

{Wqa^ = [Wq^j^ Wg^Y '^OAZ^^ 

^'^AT^ " ^^AO^^'^OAX '^OAY '^OAZ^^ " ^z^'"aT ^ 

” ^x'^'"aT " ^y*-'"AT ^ 

^'^AS^ " ^^AO^^'^OAX '^OAY Vz^^ " 

■ \^'^AS ’ V^'"as ■*■ ^AS^^®2^ ^ ^^AS^ 

" ^^SO^^'^OAX ^OAY '^OAZ^^ " ^z^^^SA^^'"aS * “aS^^V^^^x^ 

^’"SH ■" ” ^x^^^SA^^*"aS ^AS^^^y^ 

* ^’'SH " \^^^SA^^''AS ■*■ ^AS^ 
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n[rs|^ + iJ5n]te2} + {u^^} (A31) 

{Wlb> = [\o^*-'^OAX Vy ■ ^z^^\a^^*"aS ^ 

+ [T|_^][r^p][T^^][Xy] + [\f^^’"fL^*^^FA^^V^ * ^ 

■ ^y^^\A^^*"AS ‘^AS^ ^ '^\h^^’"hF^^^HA^ 

+ [TLp3[rFL][TpA] + 1 ^’'lB^^\a^^^®2^ 

- ^ . [TLH][rnp][TH,] + C\p][rpL]CTp^] 

" trLBH\A^)tv]([£g3]{u^3}) - 

+ [\^][rHp][THs][vz3C’J'y] + [\p][rp|_][Tp3][.J>^][i^y] 

+ [rLB][TLs][^z^tV)<^l> ' 

- ^^t\s^tr3H . U 3 ,] . [\H][rHp][T,3] ^ [TLp][rpL][Tp3] 

+ [rLB][TLs]Uz^)^®3> ’ + [^pK^FL^t^FA^ 

("^S ^'-s’^0 

" [^LB^[\a]) ; u- (L3.t)W [Tls]{u3,} + 

\ *^S ^ ^S ’ ^ ^ / 

+ [rLB^*^\F^^^®2^ ■ (^32) 
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From Eqs. (A27), (A32), and the fact that {Wg} = [TQ 0 ]^Wg}, we can 
write the following fully expanded expression for the inertial translational 
velocity of an arbitrary point on a main rotor blade in terms of components 
along Xj^y^Zg axes: 

" ^^DO^*-'^OAX '^OAY '^OAZ^ ’ ^z^^^DA^*-’'aS 

+ [Tgp][rpj_][Tp^][Xy][x^] + 

+ [Tgp][rp^][Tp^][Xy] + [T[)LHrLg][T|_^][Ay] 

^ ^ ' ^y^^^DA^^’^AS ^ ^AS^ 

+ t^DL^^'^LB^^\A^ ’’’ ^^DB^^*'B ■*■ ^B^^^BA^^^®2^ ■’' ^^DA^^^AS^ 

{[Tps][>^SH ‘‘' ^^DF^^*"fL^^^FA^ 

+ [TQF][rpL][Tps][i|Jz^t'f'y3 ■•■ [Tpj_][rj^g][Tj^g][t|/^][.|iy] 

+ [Tpg][rB + ^ '*' 

^ [Top][rpL][Tpg][^,] . [T^JCrLBiaLslC^,] 

+ [Tqf][*'fl][Tps] + ‘‘‘ ^^db^^’"b '*’ ‘^B^^^BS^^^®3^ 

+ [Tds]{Ush) + 3([TDp][rpL] + [TuLlCrLBlCTLp] 
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0{[Tdb 3['^B "'' ■'■ t^DB^^'^B^ 


The virtual displacement vector associated with Eq. (A17) is 

{60d} = 5XjTg^][Xy][x^]{e3} ^ ^ 

+ [Td^][^]{[^SS^^'^AS^^ ^ ^^tTBs][^^][^y]{e^} + 6^^[Tp3][^^]{e2} 

i -6v^ Hyt)] 

6u' (L5,t)V - 66[Tpp]{e2> + 6a[TpL]{e3} + 60^p[Tpg]{e^ } 
^ 6(^3 (L^jt)^ 


w’ 


+ 6<|){e^} + 6v' <{1 - v'^ - ^w'^)sin (fi^ + v'w'cos 4 >q| 
' ? 1 ? 

(1 - v' - ^ w' )cos (|)q + v'w'sin (|)q 


+ 6w' \(-l + w'^ + ^ v'^)cos + (j)sin 'I»q + ^ 4>^ cos 4 >q 

2 1 2 12 
,(-l + w' +2 V'^)sin (j)Q + (fcos ^ ^ sin <|)g 


(A34) 


Discretization of the continuous variables in Eq, (A34) via the methods 
of Sec. 3 yields 

{«0d) = «X,[Tp^][x^][x^](e3) * 

* tToA]['']([tQs^l^*AS^'^*''A* * * **y[Tps][*z](«2J 

s 


* .1, 




6e[TDF]{e2) 
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+ 6a[TpL]te3> * 


i 


w' 


+ I ■ '''^ " |-w'^)sin + v'w' 


cos 


i=l 


[(1 - v'^ - ^ w'^)cos (|>Q + v'w'sin (|iq 


'w r 

+ I Y"^ ^ \ v'^)cos (|.Q + 0 sin 'f'o + | cos 

/ 2 1? 12 
((-1 + w' + ^ v' )sin (|>Q + <|>cos 4 >q “ 2 sin «<) 



where is the matrix [$^] evaluated at the point S of the airframe. 

The first component 60^^ of { 60 ^} is needed in the main rotor virtual 
work expression of Sec. 4. It can be written from Eq. (A35) as 


*®Dx ' \zl'*^z ^ J, * %xl**x 

J ' 

$u S 

\yl% ‘ '’’dSII . 1 ^ "^Si^^S^^^i '''dS 12 ’^Si^^S^^^i 

^ ^DS13 " ^DF12 * %13 ^DBll 

^v 


(A36) 


where A^^^ is the first component of the 3 x 1 matrix [Tp^][Xy][x^]{e 2 }, 
A^^^ is the first component of the 3 x 1 matrix [Tp^][Xy]{e^ }, A^^^ Is 
the first component of the 3 x i matrix }, A^^^ is the 

first component of the 3 x l matrix [Tpg][(|;^]{e 2 }* are the first row 

components of [Tjj^][v]([jLgj][<t>^^]) , the first row, second column 

element of [Tg^]* Tqsi] is the first row, first column element of 
etc, 
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The virtual displacement vector associated with Eq. (A33) is 


{6W^} 


[TdoIC^Wq^X %AY %AZ^ ■ 


^-6v^ (Lg.t)^ 


- \ ^ fTDs]Uu3„} + {Tg}66 - {T^)6« 

6<j>S 


- ”ecp'«®cp * [To5]{6Ug) 


(A37) 


where the mat'"ices ^^Xx^’ respective coefficients of 

-X ♦ -X , etc. in Eq. (A33). 

Z X 

Discretization of the continuous variables in Eq. (A37) via the methods 
of Sec. 3 yields 

{6Wg} = t’’’pQ][5WQ^j^ 6 Wq^Y ^'^OAZ^ ■ ^''^Xz^^'^z ‘ ^''’wH 6X^} - [T,.j5X, 


+ [TqaH4»as^^^^A^ " 

S 


XX " X' 

4>x 


xy^ y 


u 

■ 1 

i=l 

4 

I 

i=l 

S. 


I Ssi'=s>«" 


i=l 



(A38) 
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t 


where is the matrix [-^j] evaliiated at point H. 

In the last term of Eq. (A38), n and c have been set equal to zero in 
antipation of the aerodynamic virtual work requirements of Sec. 4. 

The latter two components of {6Wg}, 6Wp^ and <5Wp^, are needed in the 
main rotor virtual work expression. These two quantities can be written 
from Eq. (A38) as 


6w 


Oy “ ■^D021®’'0AX * ^D022‘"0A¥ * ^D023*”cAZ ' ^z2«^z ‘ 

P P 

" jly '2j«^Aj ■ " "uS21 ,1, 

^u ^ 

' V22 '^Si^^S^'^'^ui ■ ''’uS23 ®Si .L ^2j^"’Sj 


N 

* '^62*® - \2®“ - '^ecp2®®cp ^0821 J, *gi <*)«%, ■ 
^v 

^ W2 . 1 ^ * ^0823 . 1 ^ *wi^^^%i 


(A39) 


and 


" '^0031*^'^0AX ''’d032^'^0AY ''’d 033%AZ " '’’xz3^^2 ’ '’’xx3’^\ 

P P 

■ '^Xy3'^^y ^ .1^ '^3j^"Aj ' "^3j'^’^Aj ' '^^x3^'^x ‘ \y3^'^y 

^U 

^ ^uS31 .1^ - \S32 

^ \S33 ®Si ^^S)6n^. ^ d3j.6ns. + Tg366 - T^35a 

^u \ 

* '^9cp3‘^®cp ^ "^0831 .1^ ♦ui^^^^^i ^0832 
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f 


* ^B33 ,1, *w1 


W1 


(A40) 


In Eqs. (A39) and (A40), the b^j are eleinents of the matrix [Tp^][<i>^j] , the 

elements of the matrix t^g^5K'7]([£Qg][4>^2]), and the d^j are elements 
of the matrix [Tp33C«5(^j. 

Finally, the third component of virtual displacement vector associated 
with Eq. (A29) is 

^'^ATz " ^A031%AX ^A032^'''0Ay ^A033%Az ' ^Xz3^^z “ ^Xx3‘^\ 

P 

■ " I *AT3j^’^Aj (A41) 

w * 


where is the third component of the 3 x 1 matrix 

^Xx3 component of the 3 x 1 matrix 

is the third component of the 3 x 1 matrix [r.., + u.^ 

A I AT 

are the third row elements of the matrix 


['‘at * 

* *'Xy3 

jie^}, and the 





I 
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